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Parallel Computers 


Scientists at the Naval Ocean Systems Center are ex- 
ploring the design of massively parallel computers for the 
real-time processing of sensor data. They designed an IBM- 
AT board containing 32 transputers (a British-made high- 
performance RISC CPU specifically designed for 
parallelism). A unique feature of this approach is that the 
interconnect topology between all of these processors is 
variable under software control via four crossbar-switch 
chips. Therefore, algorithms can be easily run on different 


architectures(i.e.,ring, mesh, line, hypercube, tree) toexplore 
efficiency issues. The four-board set in the photograph con- 
tains a hundred transputers, each one capable of over 10 
million instructions per second and 1.5 million floating-point 
operations per second. This combination of high perfor- 
mance, reconfigurration flexibility, and low cost is well 
suited to challenging problems in image processing, 
multi-target tracking, data fusion, radar processing, and 
underwater acoustic surveillance. 
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Introductory Remarks by 
Dr. Duston and Dr. Saalfeld 


Dr. Dwight Duston 


Director, Innovative Science and Technology 
Office of the Strategic Defense Initiative Organization 


Five years ago, the Strategic Defense Initiative Organization 
was established to spearhead one of the most ambitious research and 
development programs of the twentieth century...the search for defen- 
ses against the threat of ballistic missiles. The Innovative Science and 
Technology (IST) Office within SDIO was established soon after to 
fund novel approaches and explore new concepts in strategic defense. 
Today, with a program that sponsors over 40 basic research projects 
in a wide variety of disciplines and technologies and with an annual 
budget exceeding $150 M, the IST Program supports over 600 
principal investigators at universities, small businesses, national and 
government laboratories, and in large industry. 

The Office of Naval Research has played a key role in the growth 
and continued success of the IST Program. Through its talented corps 
of scientific officers, over half of IST’s 40 projects are technically 
managed by ONR. The Science & Technology Agents, or STA’s, have 
the responsibility for receiving and reviewing proposals, initiating 
procurement actions, managing individual research projects on a 
day-to-day basis, and providing technical and programmatic infor- 
mation on a regular basis to the IST staff members. 

While this may sound very rigorous and procedural, the relation- 
ship between IST and ONR has proven to be anything but. Both 
organizations have strived to form a strongly synergistic interaction. 
The technical opinion of the ONR STA counts heavily in the IST 
decision-making process. The relationship is more collegial than 
managerial. This unique environment that exists between agencies 
has shown to be a profitable one for both SDIO and the Navy. 
Programs leverage off of each other, stronger interaction and col- 
laboration occurs between both groups of researchers, and, in many 
cases, Navy and SDI programs have very similar goals and require- 
ments. This will be evident as you read through some of the articles 
contained herein. 

The SDIO Innovative Science and Technology Office is very 
proud of its association with the Office of Naval Research, a thriving 
cooperative effort between the Department of Defense’s oldest and 
youngest basics research organizations! As the first five years of this 
cooperation come to a close, IST looks forward to a future of 
successful and rewarding interaction between us. 
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Dr. Fred Saalfeld 


Director, Office of Naval Research 
Office of the Chief of Naval Research 


When the Innovative Science and Technology (IST) Office of 
the Strategic Defense Initiative Organization was established five 
years ago, to pursue the basic research programs of the organization, 
the Office of Naval Research (ONR) was there to lend a hand. Today 
ONR and IST are partners in basic research, benefiting IST’s mission 
to develop a non-nuclear defense against ballistic missiles, and help- 
ing ONR’s mission to prepare for the Navy of the future. 

ONR is the oldest federal agency funding basic scientific re- 
search. For 43 years it has attracted distinguished scientists, primarily 
from the universities, to participate in research programs and provide 
the basis for our future technology. Seventeen of ONR’s existing and 
long established programs were vital to IST’s major research interests 
in: high speed computing; sensing, discrimination and signal process- 
ing; space power and power conditioning; directed and kinetic energy 
concepts; materials and structures; and propulsion and propellants. 
For these programs, IST and ONR have joined forces. 

In this issue of Naval Research Reviews, we offer articles 
describing six ONR and IST’s cooperative programs. Fred Quelle in 
“Superconducting IR Focal Plane Arrays” discusses some major 
problems confronting the development of semiconducting IR focal 
plane detectors, and how new technology based on superconductivity 
may overcome them. Clifford Lau and colleagues in “Research in 
Advanced Electronic System Reliability” describes research to 
develop a scientific understanding of the factors impacting the 
lifetime reliability of microelectronic systems, including Very Large 
Scale Integration (VLSI) technology and submicron devices. Steven 
Fishman and co-authors in “Advanced Composites for Strategic 
Defense Initiative Systems” couples micromechanics, materials 
science, chemistry, and physics to explore the influence of atomic 
scale processes of composite materials on their bulk macroscopic 
behavicr. Jack Davis in “Imploding Krypton Gas Puff Plasmas as 
Potential Source for X-Ray Lasers” considers using pulsed power 
generators to drive a configuration of coaxially cylindrically sym- 
metric krypton gas puff plasmas to conditions conductive to achieving 
gain in several lasing transitions of neonlike krypton. Max Yoder in 
“Artifact Diarnond: Its Allure and Significance” discusses how alter- 
native materials, such as tetrahedrally bonded carbon films for 
microelectronic components, can enhance system performance in a 
high-stress space environment. A. Joshe describes in “The Oxygen 
Heat Engine Based on Solid Electrolytes” a promising thermoelectric 
energy conversion concept which converts heat to electricity directly 
at high efficiency. 

ONR often gets research started in critical fields without having 
to carry the total support burden. The joint and/or comanaged 
programs with IST, described here, are examples of such manage- 
ment. ONR is proud of its productive relationship with SDIO-IST. We 
look forward to continuing cooperative programs providing the finest 
scientific research for our national defense and the advancement of 
human knowledge. 





Superconducting IR 
Focal Plane Arrays 


Fred W. Quelle, Jr. 
Office of Naval Research 


Abstract 


This article discusses a number of the major problems 
confronting the development of present day semiconducting 
IR focal plane detectors for large array space based applica- 
tions, shows that an entirely new detector technology based on 
superconductivity may circumvent many of these difficulties, 
provides detailed data on the characteristics of these devices 
and outlines the development program underway to exploit 
this new technology. 


Introduction 


Present Semiconductor 
Focal Plane Arrays 


Present infrared focal plane array technology involves the 
assembly of elements consisting of a high impedance semi- 
conducting detector followed by a semiconductor amplifier, 
into a two dimensional detector array whose outputs are 
muxed together to produce an analog signal. This signal cor- 
responds to the sequential interrogation of the infrared energy 
falling on each of the individual detector elements. Usually the 
detectors are fabricated in a two dimensional array on one 
substrate and the amplifier and muxing circuitry in a cor- 
responding two dimensional array on a second substrate. Each 
of the detector and amplifier elements on these two arrays are 
then connected together by means of strategically placed 
indium bump bonds sandwiched between the two substrates. 

Even with the benefit of today’s modern large scale in- 
tegrated circuit technology, it is a nontrivial task to fabricate 
large arrays of semiconducting elements on a single substrate 
having uniform electrical characteristics with no or at worst 
only a few bad elements. The highly precise impurity compen- 
sation and doping of the semiconductor substrate that is re- 


Three/1989 3 





quired in order to construct these devices makes it difficult to 
achieve high yields and further dictates constraints on mini- 
mum component size. For this reason it is difficult to construct 
the required semiconductor amplifier and readout circuitry in 
much less than a 100 micron square. Thus for focal plane 
arrays that are appreciably larger than 100x100 elements, it is 
necessary to assemble a mosaic of detector arrays. 


System Implications Of Present Devices 


In addition to the above fabrication problems, the semi- 
conductor electronics associated with the focal plane array 
dissipates a substantial amount of power. The Strategic 
Defense Initiative (SDI) has used a figures of ~10 microwatts 
per pixel in designing their systems although currently used 
devices have significantly greater power dissipations. Even 
using the SDI figure of 10 microwatts per pixel says that a 
1000x1000 array will have a heat dissipation from its 
electronics of 10 watts. Since noise considerations require 
thermal imaging detectors to be cooled to between 40 and 100° 
Kelvin (depending on the application), this power dissipation 
produces a serious refrigeration problem. To make matters worse, 
the large number of data lines required to bring the analog signal 
information from the focal planes to the high temperature 
processing electronics, conduct significant amounts of ther- 
mal energy down to the focal plane. The drive electronics 


associated with these lines produce additional heat dissipation. 
All of this energy inust also be removed by the refrigerator. 
Figure 1 gives SDI projections of the weight and electrical 
power requirements for several types of refrigeration systems 
as a function of refrigeration power. From this, the serious 
systems implications resulting from the use of semiconducting 
IR focal plane array readout electronics should be apparent. 


Potential Of Superconducting Devices 


Since noise considerations require thermal imaging IR 
focal plane arrays to be cooled to relatively low temperatures, 
the possibility of replacing the semiconducting amplifier cir- 
cuitry with essentially dissipationless superconducting 
electronics is appealing. Such an approach could reduce the 
power dissipation of the electronics by more than two orders 
of magnitude. Even if the lower temperatures necessary for the 
superconducting electronics were to increase the weight and 
power requirements of the refrigeration equipment needed to 
provide a given amount of cooling power by as much as an 
order of magnitude, we would still be an order of magnitude 
ahead of the game because of the much lower cooling power 
requirement. Furthermore, if on chip low dissipation 
electronics were to be available, it would be possible to do A/D 
conversion and more processing on chip and come out of the 
low temperature focal plane array with a single digital output 





Figure 1 


SDI projections of the weight and electrical power requirements for several types of refrigeration systems as a function of 
refrigeration power. 
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Figure 2 
A band gap situated symmetrically about the Fermi level. 
































stream sent up an optical fiber. This would avoid much of the 
thermal conduction heating produced by large numbers of 
analog signal lines. 

Superconducting SQUID amplifiers are inherently ex- 
tremely low noise devices. They have been shown to have 
noise levels as much as two orders of magnitude below that of 
comparable semiconducting devices. It remains to be proven 
whether practical superconducting electronic circuits will ex- 
hibit this much of an improvement over the present semicon- 
ducting electronics used in focal plane arrays; however, it is 
reasonable to assume that a significant improvement in noise 
characteristics will be achieved by employing superconduct- 
ing electronics. 

Because of the inherently simple nature of monolithic 
superconducting electronics (wires and insulators placed on a 
glass or quartz substrate as opposed to carefully doped semi- 
conductors) very high production yields can be achieved. This 
technology can readily be scaled to dimensions of many inches 
and is intrinsically capable of producing large uniform arrays 
having few if any bad elements. In addition, the smaller size 
of superconducting electronic circuits (few microns as op- 
posed to many 10’s of microns for semiconductor circuits), 
permits much more sophisticated processing to be carried out 
in a given area on the focal plane. 

With all of these advantages, there is still a serious prob- 
lem. Superconducting electronics are basically low im- 
pedance devices that cannot be directly driven by high 
impedance semiconductor detectors unless impedance trans- 
formation circuitry is used. Passive impedance transformation 
would require a transformer with a high turnsratio and thus is 
not very practical for this application. Alternatively, the use of 
a semiconductor transimpedance amplifier to match the im- 
pedance would sacrifice most of the gain which we are at- 
tempting to achieve with the superconducting circuitry. What 
is needed is a low impedance optical detector whose sensitivity 
is at least competitive with present semiconducting detectors. 
Since we have already justified committing ourselves to su- 
perconducting cryogenic temperatures to support supercon- 
ducting electronics, it is reasonable to consider whether a 
practical superconducting optical detector can be developed. 


Superconducting Optical 


Detectors 


At the time the search for a superconducting optical 
detector was begun, no one had previously attempted to make 
such a detector. High temperature superconductors had not 
been invented. Furthermore, the high optical reflectivity char- 
acteristic of metallic superconductors suggested that even if 
we were able to produce a reasonably sensitive optical detec- 
tor, it’s quantum efficiency would be low due to the high 
reflectivity. Most workers in the superconductivity field were 
extremely pessimistic regarding the practicality of developing 
such a detector. 


Possible Superconducting 
Optical Detectors 


The resistivity of a simple superconductor below its tran- 
sition temperature is zero and remains there until the tempera- 
ture is raised up to the transition temperature at which time its 
resistance rapidly rises to the normal resistance of the material. 
In principal one could attempt to bolometerically use this rapid 
change in resistance to detect optical radiation. For metallic 
superconductors this transition occurs within 1 or 2 mil- 
lidegrees Kelvin so while dR/dT may be phenomenally large, 
maintaining a large array of detectors in this few millidegree 
window is impractical. Some of the new high temperature 
ceramic superconductors have a broader transition region; 
however, since their transition temperature is well above the 
transition temperature of the metallic superconductors that 
currently have to be used for the electronics, they would be 
incompatible with our present superconducting electronics. 
Thus using the superconducting to normal transition as a 
conventional bolometric detector does not appear to be a 
viable candidate for use with presently available supercon- 
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ducting electronic circuits. Some time in the future, if we 
should be successful in fabricating practical superconducting 
electronic circuits from high temperature superconductor 
materials, this issue may have to be reevaluated. 


The Physics Of Josephson 
Optical Detectors 


Superconductors have characteristics which to some ex- 
tent are analogous to a semiconductors. The supercurrent in a 
superconductor is due to the correlated motion of Cooper pairs 
located at the Fermi level in the superconductor. This corre- 
lated behavior also gives rise to a forbidden region of energy 
States (ie. a band gap) situated symmetrically about the Fermi 
level as illustrated in Figure 2. As in the case of a semiconduc- 
tor, above and below this band gap there is a very high density 
of states. In typical metallic superconductors this bandgap is 
of the order of a few milli-electron volts as opposed to 
bandgaps in semiconductors which typically are of the order 
of an electron volt. 

Josephenson devices consists of two superconductors 
which are separated by a thin insulating layer as shown in 
Figure 3. When no potential is placed across the Josephenson 
device, the two sets of supercurrent states located at the Fermi 
level in each of the two respective superconductors are just 
opposite one another and tunneling through the thin insulating 
layer gives rise to a substantial supercurrent through the 
device. When a bias potential is placed across the Josephenson 
device, the tunneling current ceases until the bias is sufficient 
to bring one of the filled bands to an energy level opposite the 
large number of empty states in the other electrode. (See 
Figure 4.) Tunneling of carriers through the thin insulating 





Figure 4 


When a bias potential is placed across the Josephson device, 
the tunneling current ceases until the bais is sufficient to bring 
one of the filled bands to an energy level opposite the large 
number of empty states in the other electrodes. 
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Tunneling of carriers through the thin insuiating layer then 
begins again giving rise to the current voltage characteristics 
illustrated here. 














layer then begins again giving rise to the current voltage 
characteristics illustrated in Figure 5. The bias potential across 
the Josephenson junction at which the sharp increase in current 
through the junction occurs is a very precise measure of the 
superconducting band gap. Thus a Josephenson device can be 
used to very accurately determine the superconducting band 
gap. 

The width of the superconducting band gap in a supercon- 
ductor below its transition temperature is related to the density 
of Cooper pairs in the superconductor both of which go to zero 
when the temperature of the superconductor is raised to the 
transition temperature. This temperature dependance of the 
bandgap is illustrated in Figure 6. Thus an optically induced 
thermal change in the temperature of a superconductor will 
result in a change of its superconducting band gap. Cooper 
pairs can also be broken through direct interactions with 
incident photons or electrons giving rise to a second 
mechanism known as pair breaking for modifying the bank 
gap of a Josephenson junction. both of these phenomenon 
could give rise to a change in the superconducting bandgap 
and thus an optical sensitivity in Josephenson devices. Pre- 
vious experience with superconductors, which in retrospect 
were thermally well coupled to the low temperature substrate 
or bath, lead us to focus primarily on the pair breaking 
phenomena in our attempts to develop an optical sensitive 
superconducting device. 


Practical Detector Embodiments 


In discussing possibilities for designing a superconduct- 
ing optical detector with various members of the supercon- 
ducting community in 1984-5, Arnold Silver of TRW was the 
first to come up with a potentially viable solution. He felt that 
the most serious problem would be the high reflectivity of the 





superconductor. He therefore proposed to absorb the infrared 
photons in a semiconductor substrate material such as GaAs. 
On the back side of the substrate he would place a Josephenson 
Junction. The absorbed photon would create an electron hole 
pair in the semiconductor and by means of a potential placed 
across the semiconductor he proposed to accelerate the 
electron and propel it into the Josephenson Junction. (See 
Figure 7.) Silver felt that the electron would have an energy 
of the order of 1000 times the bandgap energy and thus was 
bound to significantly modify the properties ‘of the 
Josephenson Junction but there was no basis on which to 
calculate the efficiency of such a process. The only thing 
one could do was to try it and see whether and how well it 
worked. TRW was given a small contract and began work 
in March of 1985. Unfortunately, problems of fabricating 
Josephenson devices on GaAs substrates while maintaining 
appropriate electrical properties in the GaAs proved to be a 
more substantial task than originally envisioned. After over- 
coming numerous technical problems data now indicates 
that current gains (change in Josephenson current per 
electron injected from the semiconductor) in excess of 140 
have been achieved. This current gain is in principal noise 
free in the same way that the current gain in a photomul- 
tiplier tube is noise free; however, detailed measurements 
to verify this have yet to be completed. This injected current 
can be gaited so that it is dumped into the Josephenson 
junction just at the time the circuit is being read out. Such 
an approach could considerably enhance sensitivity and 
signal to noise ratios in large focal plane arrays. 

In the mean time Sadig Faris of HYPRES was aware of 
Japanese work on BaPbBi which appeared to exhibit a rela- 
tively good photo response. BaPbBi did not appear to offer 
exciting possibilities for fabricating large arrays of uniform 
optical detectors, so Faris proposed to see whether he could 
achieve similar results in Nb and NDN. He felt that the high 
optical reflectivity problem could be overcome by anti-reflec- 
tive coatings on the superconductor. He was given a small 
contract in August 1985 to begin exploring these possibilities. 
The HYPRES work, which only involved laying down 
Josephenson Junctions on SiO2 substrates, proceeded more 
rapidly than the TRW program and within a year they were 
measuring responsivities (R) of 3-8000 volts/watt for 10x10 
micron junctions. Spectral response has been observed to be 
reasonably uniform from the visible out to 20 microns. 


Characteristics of Present Devices 


After demonstrating reasonable responsivity, the next 
step was to measure the noise equivalent power (NEP=Vn/R). 
For this, semiconductor amplifiers having input noise voltages 
(Vn) of 3) V/VHz have been initially used since they did not 
require the development of a squid amplifier. Amplifier noise 
in all cases dominated the output even when a 200 element 
series array of superconducting Josephenson junctions was 





Figure 6 
Temperature dependence of the bandgap. 
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used in an attempt to overcome the amplifier noise. Neverthe- 
less, NEP’s=3x10"*W/VHz have been obtained for 10x10 
micron junctions. For a 1x1 micron detectors, which are quite 
practical to fabricate, the change in bandgap with intensity will 
be the same as for the larger junction but the power intercepted 
by the detector will be two orders of magnitude lower. Thus it 
will have an NEP=3x10"°W/VHz! If pixels larger than 1x1 
micron are required, a series of 1x1 micron detectors would 
be used to fill the pixel area thereby boosting the output signal 
by N the number of elements in the series array. If when 
superconducting amplifiers are employed, the dominant noise 
turns out not to be amplifier front end noise but rather detector 
noise, this series arrangement would cause the detector noise 
to go up by only VN. 

It was initially believed that pair breaking was the source of 
the optically induced change in the superconducting bandgap of 
measured by the Josephenson Junctions. For the TRW devices in 
which the Josephenson Junction is placed on a high thermal 
conductivity GaAs substrate this is still the correct interpretation. 
The HYPRES devices have usually been fabricated on lower 
thermal conductivity glass substrates. In this case it has 
recently become apparent that the change in bandgap due to 
the optically induced heating of the Josephenson junction 
appears to be an order of magnitude greater than the contribu- 
tion due to pair breaking. When the HYPRES devices are 
placed on Si substrates or immersed in liquid helium this 
enhanced thermal effect disappears and we are left with only 
the response due to pair breaking. The thermal time constant 
for these devices on glass are in the submillisecond range and 
thus provide reasonable response times for a focal plane array. 
There are obviously many parameters that need to be ex- 
plored in order to optimally exploit this phenomena. Since 
these devices function thermally, the optical reflectivity prob- 
lem can now be solved by coating them with an appropriate 
absorbing material. A properly designed black coating should 
yield a factor of 2 or more increase in the above NEP’s. 
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Figure 7 


Accelerating the electron and propelling it into the Josephson Junction. 
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Further Development and 
Characterization Of Detectors 


The major thrust of these programs during the coming 
year will be to install low noise superconducting amplifiers 
and begin to explore the fundamental noise limits with these 
devices. To date there is little information on which to base 
speculations of what those limiting noise mechanisms will be 
or the limiting sensitivity that they will provide. Based on the 
presently available data obtained using semiconductor 
amplifiers and knowledge of the improvements in noise char- 
acteristics that have generally been associated with supercon- 
ducting circuits, it is anticipated that significant improvements 
in sensitivity will be forthcoming. 


Superconducting Focal 
Plane Arrays 


Superconducting Version Of Present 
Focal Plane Arrays 


Next year, small arrays, 4x4 and 8x8, will be fabricated 
and tested. After that the next step will be to devise appropriate 
integrated superconducting amplifying, muxing and readout 
electronics and to devise ways to monolithically integrate all 
of this circuitry on a single substrate. Hopefully, this can be 
accomplished with circuitry having a heat dissipation of the 
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order of .1 microwatt per pixel. It is not unreasonable to 
contemplate fabrication of large superconducting focal plane 
arrays having multi-inch dimensions using presently available 
superconducting fabrication techniques. The implications of 
this to IR focal plane array technology could be revolutionary. 


Development Of Additional Onchip 
Processing 


A major payoff for incorporating low heat dissipation 
superconducting electronics into focal plane arrays is that heat 
generation will no longer be the primary limiting factor on how 
much processing can be done on chip. A/D conversion is a task 
which can readily be accomplished with superconducting 
circuitry. Furthermore, digital processing associated with such 
things as target motion, etc. can be handled on chip. All of this 
digital output can be muxed into a single electrooptic 
modulator modulating a fiberoptic light beam to communicate 
with the outside world. It is anticipated that all of this can be 
accomplishable at a small fraction of the power dissipation of 
present day semiconducting focal plane arrays. 


Biography 


Dr. Fred Quelle is a scientific officer for the Applied 
Physics Field Division of the Office of Naval Research, Bos- 
ton. His publications and research have been primarily in the 
fields of lasers, electro-optics, and detectors. 





Research In Advanced 
Electronic System 


Reliability 


C. Lau, Office Of Naval Research 
C. Hu, University Of California, Berkeley 
E. J. McCluskey, Stanford University 


Abstract 


Reliability of microelectronic chips and systems is be- 
coming an increasingly important factor in the quality of 
electronic systems as the scale of integration continues to 
increase while device geometry continues to shrink to smaller 
dimensions. In the case of Strategic Defense Initiative (SDI) 
systems, designed-in reliability is imperative because of the 
criticality of the electronic systems and the impossibility of 
maintenance and repair in space. Microelectronic systems for 
sensors, signal processors, and supercomputers must be able 
to operate without human intervention for an extended period 
of time in a very harsh space radiation environment. Four 
specific areas of research in the reliability/survivability of very 
large scale integration (VLSI) circuits have been identified as 
Critical to the success of future microelectronic systems for 
SDI application. These are fault tolerant systems and architec- 
tures, thin dielectric insulator films, contact/interconnect 
metal electromigration, and radiation hardening. Research has 
increased our understanding of the physics of dielectric break- 
downs and wearouts and the mechanism of charge trapping in 
thin gate oxides and the understanding of the electromigration 
process in thin-film metals. Research in system reliability will 
continue to develop multiple fault models and methods for 
self-testing, error detection, and error control. Because 
reliability is a system-level issue, hierarchical methodology 
must be developed that spans reliability research in solid state 
materials, devices, circuits, macrocells, chips, packaging, 
boards, subsystems, and systems. 


Introduction 


Although not often explicitly stated, reliability has always 
been a major concern of electronic system designers. 
Reliability was very important in the days of vacuum tubes. 
Because of the deterioration of the heating filament and the 
electron emitting cathode, a vacuum tube typically changed its 
characteristics such as transconductance over a period of time, 
and lasted only for a few years. In fact, the very first digital 
computer built with vacuum tubes and mechanical relays did 
not become operational until the tubes could function reliably. 
Thus circuit designers had to design systems with a large set 
of varying devices, and had to consider reliability in the 
design. 

Since the invention of the transistor, the reliability picture 
began to change. The transistor, once perfected, proved to be 
extremely reliable. Since there was no heating element or 
cathode, the transistor could virtually last forever, except for 
the gradual wearout of the solid state material. Thus systems 
built with transistors and solid state devices could last for over 
a decade without deterioration in performance. 

The integrated circuit (IC) was made by putting many 
transistors monolithicly on the same substrate. For small scale 
integration circuits (SSI), with hundreds of devices on a chip, 
the reliability was not too different from discrete transistors 
except for the fact that there was greater power dissipation and 
consequently greater temperature rise on the chip. In other 
words, SSI integrated circuits, once tested to be working 
properly, could last for years without failure. Meanwhile the 
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scale of integration continued to increase and the device 
geometry shrank to ever smaller dimensions. Today it is not 
uncommon to see very large scale integration (VLSI) circuits 
containing a million devices at one micron geometry. There is 
every reason to believe that this trend to increase the scale of 
integration will continue. 

As device geometry shrinks to submicron size and the scale 
of integration goes beyond VLSI, the reliability picture again 
begins to change. The once reliable discrete transistor is now so 
small that any slight change in the characteristic of the silicon 
material would cause it to fail. Furthermore, in digital circuits 
with hundreds of thousands of gates, the failure of a single device 
would cause the entire VLSI chip to fail. Thus as microelectronics 
become more and more complex, there are serious concerns 
about the reliability of the VLSI chips and systems. There are 
several factors which contribute to these concerns. 


Scaling 


All insulating materials breakdown if given a sufficiently 
high electric field. So it is true with the gate dielectric in 
Complementary Metal Oxide Semiconductor (CMOS) tech- 
nology. All gate dielectrics in today’s CMOS VLSI circuits are 
silicon oxides (SiO2), which break down at about 10 MV/cm 
in a few seconds. Breakdown occurs at a lower electric field 
if the field is applied for a longer period of time and if latent 
defects are present. As the device geometry shrinks to below 
a micron, the thickness of the gate dielectric also shrinks to 
about a hundred angstroms. If the supply voltage remains at 5 
volts, the electric field at the gate is 5 MV/cm, which is near 
or above the practical breakdown strength of SiO2. The electric 
field along the electron path in the silicon near the drain also 
becomes quite high. This high field can create “hot electrons”, 
which can damage the Si/SiO> interface and cause the transis- 
tor performance to degrade with time. These are the primary 
reasons for the push to lower the supply voltage to about 3 
volts for submicron circuits. 


Complexity 


As the scale of integration continues to increase, more and 
more devices are being put on a chip. One reason is that it is 
more efficient utilization of silicon to put related circuits on 
the same chip. Another reason is that circuits on the same chip 
can run at a higher clock speed than off-chip circuits because 
of shorter interconnect distances with lower resistances and 
capacitances. Off-chip driver circuits are typically large and 
consume more power, and are to be avoided. Unfortunately in 
digital circuits without any designed-in fault tolerance, if one 
of the devices fails the entire chip is useless. With that many 
devices on a chip, the possibility that one of them will fail is 
not minuscule. 
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Speed 


There is a continuing demand for faster and faster chips. 
For civilian application, the clock speed on many 
microprocessor chips have been pushed from a few megahertz 
(MHz) to as high as 20 MHz in the last few years. The clock 
speed on the DoD Very High Speed Integrated Circuits 
(VHSIC) phase 1 program is 25 MHz for 1.25 micron technol- 
ogy. For the VHSIC phase 2 program, the on-chip clock is to 
be 100 MHz for the 0.5 micron technology. But as the clock 
frequency increases, there is greater susceptibility to coupling 
between wires on the chip. At high frequencies, the short wires 
would behave like antennas and radiate out electromagnetic 
(EM) energy. This EM radiation gets picked up by the wires 
in its vicinity. This coupling contributes to greater noise on the 
chip. At the same time as the devices become smaller, there is 
a need for greater precision in the threshold voltages. Thus the 
noise margin, or the ability of the device to withstand noisy 
signal without falsely switching states, is lowered and causes 
the chip to be susceptible to errors. 


Process Variations 


As the devices become smaller and the chips get bigger, 
the fabrication and manufacturing process is also getting more 
complex. Increasingly exotic devices are invented in that 
endless pursuit for performance and reliability. For example, 
lightly doped drain (LDD) devices are widely used to alleviate 
the concentration of electric fields near the drain areas of Field 
Effect Transistors (FET), and vertical trench isolation techni- 
ques have been tried to prevent the coupling of closely located 
devices. Multilevel metalizations are also used to more effi- 
ciently utilize the silicon areas on the chip and to reduce the 
overall size of the chip. But multilevel metals require more 
masking layers and thus further increase process complexity. 
As the devices are so small, the process parameters must be 
strictly controlled to maintain the uniformity of the fabricated 
devices. Any deviations from the designed parameters and 
geometry can result in a bad chip. For example for 1 micron 
geometry, if the fabricated MOSFET channel length is 0.2 
micron longer than the design value, the circuit will probably 
fail the speed specification. On the other hand, if the channel 
length is 0.2 micron shorter, the circuit may fail after a short 
service time due to hot carrier induced degradation. Further- 
more, when the FET channel length is below a micron, the 
smallest particle can contaminate the process and cause the 
fabricated device to fail. Thus the control of purity of the 
chemicals, gases, water, and clean-room air is becoming in- 
creasingly critical to the yield and reliability of the fabricated 
VLSI chips. 





Material Defects 


All bulk materials contain defects regardless of the great 
care that is put into its preparation. As devices get smaller, the 
smallest defect in the silicon material or metalization can cause 
malfunction in the devices and circuits. At the same time, as 
the VLSI chips become bigger, it becomes critical to reduce 
the number of defects per area since the yield is exponentially 
lowered as a function of the chip size for a given defect density. 
However, in spite of significant research efforts, it is extremely 
difficult to grow silicon crystals without any defects. The best 
silicon substrates still may contain as many as five defects per 
centimeter square. This would result in very low yields (<1%) 
for extremely large chips (1 cm”). Furthermore, even if the chip 
is tested to have functioned properly initially, latent defects 
may cause reliability problems later on in the chip’s lifetime. 


SDI/IST Program 


Reliability is important to military electronic systems. 
The Navy ships and aircrafts must be able to sense the threat 
environment, develop defense strategies, and utilize its 
weapons to negate that threat. The weapon systems and its 
electronics must work when it is needed in the battle environ- 
ment. There is no second chance. Compounding the reliability 
problem is that the military environment is much more severe 
in terms of temperature ranges, shock and vibration, humidity, 
and most importantly radiation. In fact the radiation require- 
ments, such as total dose, dose rate, and single event upsets are 
getting harder and harder to meet. 

For strategic defense systems, the reliability of VLSI 
chips and systems is particularly important. The importance of 
reliability of electronic systems to the Strategic Defense Initia- 
tive (SDI) is outlined in the Fletcher reports.'?_ VLSI chips 
and systems are undoubtedly the underlying technology that 
makes large-scale strategic defense system possible, including 
sensors, signal processors, and supercomputers for battle 
management. Because of the criticality of SDI systems, the 
electronics must operate continuously and particularly at time 
of high demand. The lifetime of space-based systems must be 
ten to twenty or more years. Unfortunately, maintenance and 
repair are virtually impossible for space-based electronic sys- 
tems. Thus designed-in reliability is imperative for these very 
complex electronic systems. Additionally, the massive nature 
of SDI systems (thousands of sensors and parallel computers) 
would require generic reliability techniques to make it affor- 
dable. It is with these considerations that the Innovative 
Science and Technology (IST) office of the SDI Office in- 
itiated a research program on Reliable/Survivable Advanced 
Electronic Systems in 1985. 

The objective of this research program is to develop the 
scientific understanding of factors that impact the lifetime 





Figure 1 
Hierarchial Reliability 





MATERIALS 





reliability of microelectronic systems, and to develop the 
methodology for the design of highly reliabie/survivable ad- 
vanced electronic systems. The program is in direct support of 
virtually every aspect of the SDI program, particularly in the 
areas of sensors, signal processors, and battle management. 
From the very start, it is recognized that VLSI reliability is a 
very complex problem. The approach is to develop a 
methodology of hierarchial reliability, spanning from the 
purity of starting silicon materials to the design of reliable 
systems, as shown in Figure 1. The methodology must ensure 
that errors that are not detected at a lower level should be 
detected at a higher level. For example, error detection and 
correction circuits must be included at the microcell level to 
correct for faults that may go undetected at the circuit level or 
the device level. Similarly self-testing and reconfiguration 
schemes should be developed at the system level to ensure 
total system reliability. However due to funding limitations, it 
is not possible to fund all areas of research. Four areas of VLSI 
reliability have been identified as critical to the success of 
future VLSI chips and systems for SDI application: (1) Sys- 
tems and architectures for fault tolerance; (2) Thin dielectric 
insulator films; (3) Contact/interconnect metal electromigra- 
tion; and (4) Radiation effects on VLSI devices and circuits, 
and the associated radiation hardening techniques. 
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Technology Reliability Issues 


Due to the limited space we will only discuss the technol- 
ogy reliability issues relating to metal-oxide-semiconductor 
(MOS) technology — the workhorse of VLSI. The need for 
increased circuit density and speed has driven down the lateral 
dimension of the transistor (the channel length ) and the 
vertical dimension of the transistor (the gate oxide thickness). 
This scaling results in larger transistor currents, which are 
responsible for the improved circuit speed. 


Oxide 
thickness 
1200A 
700A 
250A 
200A 
150A in DRAM 
80A in EEPROM 


Current 
per mm 


15mA 
50mA 
200mA 
400mA 


Channel Operating 
length voltage 


7pm 5V 
4um 5V 
2 um 5V 
1 um 5V 


After such aggressive reductions in oxide thickness and 
channel length, without a corresponding decrease in the 
operating voltage, and such a large increase in current density 
itis no surprise that we are now faced with three serious threats 
to IC reliability: oxide integrity, channel hot-electron induced 
device degradation, and metal electromigration. A fourth 
reliability concern is the radiation effects. 

For each reliability problem, research must be done to 
address the underlying physical mechanisms, testing methods, 
process improvements, and to improve the reliability predic- 
tion at the circuit design stage. In the next several sections we 
shall describe each of the four problems, examples of research 
accomplishments, and research that needs to be done. 


Oxide Reliability 


Although numerous models for breakdown have been 
presented, there has not been a strong consensus concerning 
the breakdown mechanism. Electric field or current enhance- 
mentat some weak spots due to sodium ion drift, or hole trapping, 
or electron trapping, or trap generation have been suggested as 
the causes of breakdown at various times in the past. 

Recent research has uncovered strong evidences of hole 
generation and trapping in the oxide during high field stressing 
and the rates of hole generation and trapping are directly 
correlated with oxide lifetime*. Hole trapping in the oxide 
Causes an increase in the electron tunneling (Fowler- 
Nordheim tunneling) current density, further increasing the 
rate of hole generation and trapping. This positive feedback 
loop leads to catastrophic breakdown‘. 
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The lifetime of a given oxide sample is a function of the 
applied electric field, a phenomenon known as time-depend- 
ent dielectric breakdown (TDDB). This makes it difficult to 
test the oxide reliability and determine a safe operating voltage 
for a given IC technology. Up to now, the standard practice is 
to assume 


tap=Ae PE (1) 


The time-to-breakdown, tgp, is measured at a few large 
oxide field, Eox, values and the test data are extrapolated 
according to Eq. (1) to determine the tgp for the intended lower 
operating field. 

The hole-induced-breakdown model described above has 
recently predicted a different relationship‘. 


tpp=Be *™ (2) 





Oxide lifetime before breakdown occurs, tgp, is a function of 
the oxide field, Eox. Lifetime at a low operating field is usually 
extrapolated from high field test data according to 
tap a &°§* However, tap 0. evo as predicted by a recent 
theory, is more accurate. 
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A method to determine the density of defects (thin spots) in 
oxide as a function of the defect severity, in terms of the 
“effective” thinning, was recently developed. This information 
allows one to predict the statistical rate of failure due to oxide 
breakdown. 
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Figure 2 shows that Eq. (2) is amuch more accurate model 
than Eq. (1). 

Oxide failure in a large VLSI circuit always occurs at a 
single weakest spot. And the reliability of the circuit is deter- 
mined by that single, randomly distributed “defect”. It is, 
therefore, imperative to characterize these “defects”. Recent 
research has developed a method to determine the defect 
density D(AX,,) as shown in Figure 3, where AXo, is the 
effective thinning of the oxide at the defect points’. Once the 
defect density for a given oxide technology has been deter- 
mined, the failure rate of a circuit containing any area of oxide 
can be predicted as a function of voltage and time according 
to Eq. (2). 

A large portion of the defects is believed to be a direct 
result of the metallic precipitates in the silicon substrate that 
become incorporated into the oxide during thermal oxidation. 
Good pre-oxidation gettering has been shown to reduce the 
rate of early failures. A novel and interesting idea, recently 
verified experimentally, is that good CVD (Chemical Vapor 
Deposition) thin oxide may be more reliable than thermal gate 
oxide because the gaseous sources contain less metallic im- 
purities than the silicon substrate. 

These university research results need to be verified 
against high volume reliability testing typically performed in 
the industry. Temperature dependence needs to be added to the 
model. The mechanism for hole generation in oxide is not 
understood, especially for the case of low voltage operation. 





Figure 4 


Inset shows the LDD MOSFET structure, which will be standard for 5V submicron or micron MOS technologies. Peak channel 
electric field (a function of the n region doping concentration) is reduced in an LDD device and the hot-carrier limited device 


lifetime can be greatly prolonged. 
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A close correlation between oxide reliability and radiation 
hardness is suspected, because both are linked to hole trapping, 
and should be investigated. 


Hot Carrier Induced MOSFET 


Degradation 


Near the drain junction in the MOSFET channel, there is 
a high field region, about 2000A long, known as the pinch-off 
region. Field strength can reach 3X10°V/cm in 1.5 um 
MOSFETs. At this electric field, electron temperature is 
2500°K. High energy electrons can enter the oxide and become 
trapped or create acceptor-like interface traps and introduce 
interface trapped charge. The result is a decrease in device 
current and a drift in the threshold voltage and other device 
parameters’. This is known as hot carrier (electron) induced 
degradation. It is customary to consider the transistor to have 
failed when the current has decreased by 3% or 10%, or the 
threshold voltage has drifted by 10 mV or 100 mV. 

The electric field strength can be significantly reduced, 
hence the device lifetime dramatically improved, by using the 
lightly doped drain (LDD) structure. An LDD structure is 
shown in the inset of Figure 4. Next to the heavily doped drain 
(and source) diffusion, a lightly doped region is introduced. 
This reduces the peak electric field in the same manner a 
decrease in doping concentration would reduce the peak field 
at a reverse biased p-n diode junction. Figure 4 shows the 
reduction in the electric field as predicted by an analytical 
model’, 

Testing and predicting device lifetime is greatly 
facilitated by noting that the hot electrons not only cause 
device degradation but also causes impact ionization. Impact 
ionization produces a small but easily measurable hole current 
that flows out of the substrate contact, known as the substrate 
current, Isup. It was theoretically predicted and now well 
established that 


—3 
hen (3) 
t=10° ~ ) 
where T is the device lifetime in seconds; W is the device width 
in 1; and Isup is in LA. This relationship is shown in Figure 5. 
Thus, t, which may be longer than 10 years, can be estimated 
by making a single I,.» measurement or a few short accelerated 
lifetime tests. 

Further research is needed to determine if device lifetime 
can be related to Isup(t) under switching voltage waveforms. If 
it can be related to Isup(t), we would be able to simulate Isup(t) 
at the circuit design stage and predict the device lifetime. If 
(and there are some evidences for this) it cannot, the reason 
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Hot-carrier limited device lifetime, t, is a power-law function of 
the impact ionization generated substrate current, Isub. W is 
the transistor width. At high voltage, for example, Isub is large 
and % is short. Such a plot allows one to extrapolate high 
voltage test data to determine the device lifetime at a lower 
operating voltage. 
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needs to be understood and alternative means for lifetime 
prediction must be found. Another challenging research area 
is how to make the interface more resistant to hot carriers. 
Hydrogen-free oxide and oxynitride are good candidates for 
hot-carrier resistant dielectrics. 





Metal Electromigration 


Electrons can knock metal (principally aluminum in VLSI 
circuits) atoms out of their lattice sites. At very high current 
density (>10°A/cm?) and over an extended period of time, 
voids and breaks in the metal line can develop. This is the basic 
electromigration problem. The mechanism is fairly well un- 
derstood and the following mean-time-to-failure model is 
universally accepted for DC stressing 


MTF = AJexp (E,/kT) (4) 


where n is about equal to 2. 

The remaining question is about pulse current stressing. 
In VLSI, current waveforms are highly peaked (CMOS in- 
verters conduct current only during the short switching tran- 
sients) at high clock frequencies (10-100 MHz). Should J in 
Eq. (4) be the average or the RMS current? Experimental 
results suggest that it is neither; MTF appears to be a compli- 
cated function of frequency and duty cycle. The effects of 
non-rectangular current pulse shapes are too intimidating to 
have been studied. 

This complex behavior is recently explained’ by a simple 
model. The metal vacancy density, n(t) is determined by the 
rate of vacancy generation, that is proportional to J(t), minus 
the rate of annihilation that is n/t, where t is the vacancy 
lifetime. The metal mass transport rate is proportional to the 
see J(t)n(t). Therefore, MTF is the time for the integral 


J-ndt to reach a certain constant. 


Future research will have to address the question of bi- 
polarity current stressing, i.e. to what extent does the current 
flowing in one direction repair the damage done by the current 
flowing in the other direction. There is also a large number of 
metal technology issues to be studied — copper doped 
aluminum, refractory metal, metal/semiconductor interface 
and metal/metal interface reliabilities. 


Radiation Effects 


There are three important radiation effects — total dose 
effect, dose rate effect, and single event upset (SEU). The total 
dose effect results from hole trapping in the oxide which 
causes the threshold voltage to drift until the circuit ceases to 
function to specification. Radiation hard processes have been 
developed to minimize this effect, but further research is 
needed. For example, some radiation hard processes attain the 
radiation hardness (unknowingly perhaps) by compensating 
the trapped holes with negative interface trap charges. Under 
certain conditions, such as at high temperature where holes can 
be thermally detrapped, this charge balance can be lost. Better 
understanding of the mechanisms and tempral dependence of 
electron and hole trapping, interface trap generation, electron 


and hole recombination, and interrelationships among these 
mechanisms is needed for a complete understanding of the 
total dose hardness. This significant overlap among radiation 
effects, hot electron effects, and oxide TDDB makes an unified 
study of these effects particularly useful and promising. New 
materials such as CVD oxide and new processes need to be 
investigated. 

The basic mechanisms of SEU and the dose rate effect are 
known and somewhat similar. Simple models have been 
proposed for new phenomena such as funneling as due to ot-par- 
ticle or heavy ions have been proposed!°. Charge collection may 
be minimized by using clever doping profiles, perhaps produced 
with high energy (MeV) ion implantation. CMOS offers addi- 
tional opportunities to manipulate the device structure. To guide 
these development, better physical understanding and analytical 
models of charge collection in complex device structures are 
needed. Circuit configurations can be modified to reduce sen- 
Sitivity to collected currents and charges. Several new SRAM 
cells have been proposed for this purpose. For designing more 
complex circuits, circuit simulators that can simulate the effects 
of, and predict the resistance to, pulse radiation and SEU would 
be extremely beneficial. 





Table 1 


Types of Tests 





Tests of manufactured parts to sort out those 

that are faulty 

Test of each die while still on the wafer. 

Test of packaged chips and separation into 
Classes or bins (military, commercial, 
industrial) 

A test to demonstrate the degree of 
compliance of a device with purchaser's 
requirements. 

Test of some but not all parts 


Test to determine whether device meets 

spcifications 

Test to determine actual values of device AC 
and DC parameters and the interaction of 
parameters. 

Used to set final specifications and to identify 
areas to improve process to increase yield. 

Test with stress (high temperature, temperature 
cycling, vibration...) applied to eliminate 
short life parts 

Reliability (Accelerated Life) Test after subjecting the part to extended high 

temperature to estimate time to failure in 

normal operation 


Test to locate failure site on failed part 


Test by quality assurance department of a 
sample of each lot of manufactured parts. 
More stringent than final test 

On-line testing to detect errors that occur 
during normal system operation 

Verifying the correctness of a design 


Productjon 


Wafer Sort or Probe 
Final or Package 


Acceptance 


Sample 
Go/No Go 


Characterization or 
Engineering 


Stress Screening (Burn-in) 


Diagnostic (Repair) 
Quality 


On-line or checking 


Design Verification 
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Figure 6 
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System Reliability Issues 


Manufacture of an electronic system that can be guaran- 
teed to operate correctly for a lifetime measured in years 
requires that: 

1. The components from which the system is assembled 
are defect-free; 

2. The assembly process is correct; 

3. The reliability of the components is sufficiently high so 
that the number that fail during the system lifetime can be 
accurately bounded; and 

4. The system includes facilities to detect the occurrence 
of a failure, isolate its effects, and reconfigure the system to 
exclude the failed device and replace its function. 

In very general terms, two complementary techniques 
must be used to ensure that these requirements are satisfied: 
the manufacturing, assembly, and design must be as fault-free 
as possible, and adequate test procedures must be used to 
detect failures. Both techniques must be used since neither is 
sufficient by itself. Producing reliable components was dis- 
cussed in the preceding section. This section discusses testing. 


Manufacturing Test 


The most important requirement in building a reliable 
system is to ensure that it is constructed of non-failed com- 
ponents. To do this, the manufacturing process should have the 
highest possible yield (ratio of good components to all com- 
ponents manufactured), and then should be followed by a 
thorough test procedure to eliminate any bad devices. 
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The current state-of-the-art in integrated circuit manufac- 
turing is that the yield can be quite low (less than 50%) for 
aggressive designs using large die with many small-geometry 
devices. Even more conventional designs rarely have yields as 
high as 90%. A board assembled with 100 chips, each having 
a 10% probability of being bad, will have a probability of 
containing a failed chip of 1 -(.9)'® = 99.997%. Clearly this 
is unacceptable; the manufactured devices must be tested to 
increase the quality level (percentage of good chips) of the 
chips used in board assembly. 

Present integrated circuit production test techniques have 
two major shoricomings: 


1. One is the use of the single-stuck fault model to rate the 
effectiveness of a test procedure, even though it is clear that 
many of the defects that actually occur are not accurately 
represented by such a model; and 

2. The other is the acceptance of a fault coverage (per- 
centage of faults detected) less than 100%. 


Fault coverage is discussed in the following subsection, 
followed by a discussion of fault models in the next subsections. 

Fault Coverage: After a piece of equipment is manufac- 
tured or assembled, it may require testing before beirg sold 
or used. The need to test depends on: (1) the process yield, Y, 
the proportion of finished units that are not defective, and (2) 
the acceptable quality level, AQL, the planned minimum 
fraction of defective shipped units. 

If, as is often the case, the process yield is lower than the 
AQL, defective units must be identified and removed so as to 
increase the percentage of good units (Figure 6). The proce- 
dure of identifying defective units is called testing (the process 





of determining whether a piece of equipment is functioning 
correctly or is broken or faulty). 

Testing done during manufacturing is called manufactur- 
ing or production testing. Testing is also done for other pur- 
poses or at other stages in a product’s life. Table 1 lists some 
of these other types of tests.!!-!2 

Test quality is usually measured by fault coverage C, the 
percentage of possible defects that are detected by the test 
procedure. The process yield and the AQL determine the 
required fault coverage. This dependence can be quantified by 
assuming a model for how manufacturing defects affect the 
process yield. The model used here!* '4 is based on the two 
papers by Williams. It assumes that the causes of chip failure 
can be enumerated as n possible defects per chip, and that the 
test process can detect the presence of m out of all n possible 
defects. It is further assumed that the occurrence of each of the 
n defects is independent and equally likely with probability p. 
The quality level achieved by a test with fault coverage C for 
a process of yield Y is given by: 


QL=y"* (5) 


This expression makes no assumptions about the range of 
values of the variables. Modern integrated circuit manufactur- 
ing typically requires quality levels of at least 99.99% (cor- 
responding to fewer than 100 defects per million devices). If 
a high yielding process (Y > 90%) is assumed (since the 
required coverage will be higher for a poorer process), it is 
possible to rewrite the expression for QL as:!° 


1-QL=(1-C)(1- Y) (6) 


For a QL of 99.99% and a 95% yield, the required fault 
coverage, C, is 99.8%. 

There is amazingly little empirical data on the relationship 
between fault coverage and quality level for integrated cir- 
cuits. Only one experiment to relate test coverage and quality 
level in an actual IC process has been reported in the litera- 
ture’®, The chip designer’s test patterns had a 96.6% fault 
coverage, but this test failed to discover 103 bad dies. Use of 
this test with 96.6% fault coverage would have resulted in an 
unacceptably low QL of 99.15% in manufactured ICs. Both 
the theory and the experiment lead to the conclusion that 
96.6% fault coverage is inadequate. 

In summary, the theoretical studies, as well as the one 
experiment that has been published, indicate that single-stuck 
production test fault coverages of greater than 99% are 
desirable for manufacturing high quality ICs. 

It is both difficult and expensive to guarantee such high 
fault coverage. Design-for-testability techniques are used in- 
creasingly to enhance the fault coverage. But even with the 
inclusion of test scan paths high fault coverage is expensive. 


One technique that has been used successfully is to partition 
the chip and test combinational subcircuits exhaustively 
(apply all possible input patterns)'’. This technique is very 
effective, but its use has been limited because of the difficulty 
of obtaining a viable partitioning. Research is on-going to 
develop automatic methods of partitioning designs.'* 

Even 100% single-stuck fault coverage may not be ade- 
quate to guarantee defect-free chips. Many failure modes are 
not correctly represented as single-stuck faults. Other fault 
models are discussed in the next subsections. 

Complex Stuck Faults: Single-stuck faults may ac- 
curately represent physical failures in which a single circuit 
net is either open-circuited or is shorted to a power rail. If, as 
is often the case, more than one such fault occurs; a multiple- 
stuck fault model is required for accurate modelling. A bridg- 
ing fault model is needed for those very common situations in 
which two signal nets are incorrectly interconnected. There is no 
disagreement about the occurrence of failures that require com- 
plex stuck fault models, multiple-stuck faults or bridging faults, 
to represent them accurately; but these models are never used in 
developing test patterns because there are so many such faults. In 
a circuit with n nets, there are approximately 2n single-stuck 
faults, 3" multiple-stuck faults, and n* two-net bridging faults. 

Both theoretical studies and experimental results have 
shown that single-stuck fault test sets detect many complex 
stuck faults. But recent studies have demonstrated that even 
test sets with 100% single-stuck fault coverage may have 
inadequate coverage of complex stuck faults'®”°. Cost-ef- 
fective techniques have been developed to augment single- 
stuck fault test sets to improve the complex fault 
coverage!*”°, Modified exhaustive test techniques, called 
pseudo-exhaustive tests, have been described that provide 
guaranteed coverage of both single and complex stuck 
faults without the necessity of fault simulation to determine 
coverage". 

Pattern Dependent Faults: There are other failure 
modes that cannot be represented as stuck faults since they do 
not result in incorrect interconnections, but rather introduce 
spurious circuit elements”. These failures are particularly 
pernicious since the element parameter values may depend 
on the chip temperature or supply voltage value and thus 
may appear to be temporary malfunctions. While stuck- 
fault test sets may fail to detect such failures, pseudo-ex- 
haustive tests have been shown to be effective in 
discovering their presence. Research is continuing to quan- 
tize the coverage of such faults. 

Delay Faults: Failures that alter component values 
without introducing extra components can cause the circuit to 
fail due to violation of timing parameters. Such failures are 
modelled as delay faults. They can be very difficult to detect 
since the parameter values can depend on temperature or 
supply voltage, and many signal paths must be tested to 
discover the presence of a localized delay fault. The develop- 
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ment of cost-effective test methods for delay faults is an area 
of intensive current research. 

MOS Faults: The possibility of having failures in NMOS 
or CMOS circuits that convert a combinational circuit into a 
sequential circuit has been demonstrated and a great deal of 
effort has been devoted to developing test procedures for such 
faults. None of these procedures are cost-effective. Further- 
more, it has been demonstrated that this is a class of failure 
that is better avoided by proper design (layout) methods”. 


Fault Tolerance 


The most critical feature of systems that must continue to 
operate correctly after a failure has occurred is the ability to 
detect and localize the failed element. Many schemes for 
error-detection have been developed, but they are all based on the 
assumption that failures manifest themselves as either single- 
stuck faults or as unidirectional errors. There is general agreement 
that most operational failures are the result of temporary rather 
than permanent component failures. Recent experiments have 
demonstrated that temporary failures do not have the charac- 
teristics of stuck faults, but rather are delay or pattern dependent 
faults . Research is underway to characterize the ability of 
standard error detector structures to detect such failures”. 


Summary 


Since the initiation of the research program and par- 
ticularly during this past year, significant accomplishments 
were made in terms of understanding the physics of failures in 
the devices and metals and developing design techniques to 
circumvent these failures. In this paper, we presented briefly 
the issues and some recent research results in technology 
reliability and system reliability. In other funded efforts, thin 
aluminum and aluminum-alloy films have been prepared and 
characterized by means of TEM, SEM, and optical micros- 
copy. Low temperature deposition of multilayer oxynitride 
thin film has been demonstrated by means of remote enhanced 
plasma chemical vapor deposition (REPCVD) technique. 
Parallel test techniques have been developed for testing high 
density memories. Reconfiguration schemes have been 
developed to ensure graceful degradation and fault tolerance 
in processor arrays. High speed high reliability signal process- 
ing using logarithmic and residue number systems have been 
demonstrated in a new experimental computer architecture. 

There are four major technology reliability concerns: 
oxide time dependent breakdown, hot carrier induced MOS- 
FET degradation, metal electromigration, and radiation ef- 
fects. For each of these four concerns we need to research these 
inter-linked subtopics — underlying physical mechanisms, 
testing methods, process improvement, and means of predict- 
ing reliability at the circuit design stage. At the present time, 
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the basic underlying mechanisms are known and rudimentary 
testing methodologies, at least for DC stressing, for comparing 
technologies are available. How to analyze tests under realistic 
waveforms and how to use the data of simple tests to predict 
circuit reliability are the next research challenges. Computer 
aided design (CAD) tools will have to be developed to ac- 
complish these research goals. There is still much promise in 
the search for improved processing technologies, material, 
and/or device structures to improve electronic component 
reliability by orders of magnitude. 

In order to achieve high system reliability and availability 
it is necessary to screen out failed parts before assembling the 
system. This requirement mandates the use of extremely 
thorough IC test techniques that detect 100% of the failures 
that can be modelled as single-stuck faults as well as close to 
100% as possible of failures requiring more complex fault 
models. This level of coverage will require the use of pseudo- 
exhaustive test techniques. CAD implementations of such 
techniques are currently being developed at the Stanford Cen- 
ter for Reliable Computing (CRC). 

Continual system operation depends on the ability to 
detect failures during operation and to reconfigure the system. 
Present techniques for failure detection are of questionable value 
for non-single-stuck faults caused by temporary failures. Work is 
underway at CRC to develop both more accurate models and also 
improved error detectors for temporary failures. 

It is clear that reliability is a very complex problem. For 
large-scale VLSI systems, reliability must start from the basic 
starting silicon materials, and going all the way to reliability 
at the system level. Since reliability is a system level issue, 
hierarchical reliability methods must be developed to solve the 
reliability problem. Reliability issues must be addressed 
throughout the manufacturing and operation of VLSI chips 
and systems at various levels, from silicon materials, fabrica- 
tion process, devices, logic design, circuits, macrocell, chips, 
boards, subsystems, and systems. The reliability problem 
cannot be solved at any single level. For example, errors that 
are not detected and corrected at the gate level must be 
detected at the macrocell level. Any remaining errors must be 
detected and corrected at the system level. 

The solution to the reliability problem lies in reliable 
devices, correct design, and thorough testing. Firstly, the sys- 
tem must be built with reliable parts since it is impossible to 
build a reliable system out of unreliable components. That is 
the reason so much research is concentrated on improving the 
reliability at the device and component level. Secondly, the 
design must be correct to make it possible to separate the 
fabrication flaws from potential design errors. Therefore, 
design techniques and CAD tools must be developed to insure 
the correctness of design before committing to silicon fabrica- 
tion. Finally, the chip and system must include designed-in 
testability to allow for the thorough testing of the system to 
insure reliable operation throughout the system’s lifetime. 
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With Office of Naval Research support, Dr. Faris created 
the first manufaciuwting facility for the production of NbN 
integrated circuit chips. These chips will be used for the 
Navy’s Digital Magnetometer project for ASW and other 
activities and HYPRES’ own products which utilize inexpen- 
sive, reliable and compact closed-cycle refrigerators operating 
below 10°K. 

Also with ONR support, Dr. Faris, demonstrated for the 
first time that superconducting sensors have higher sen- 
Sitivities than conventional technologies. The supercon- 
ducting sensors will be the major component of monolithic 
focal plane arrays which combine the sensor and other 
low-power dissipation, superconducting circuits for digital 
and analog processing. The focal plane arrays will exhibit 
higher performance, lower cost, more reliability and lower 
weight. 
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Future systems for the Strategic Defense Initiative have 
material requirements which are more severe than for other 
Department of Defense or industrial applications. These re- 
quirements include operation at extremely high temperatures, 
absolute dimensional stability during thermal cycling, high 
vibration damping, low weight, high strength and mechanical 
stiffness, and high reliability. Added to these is the necessity 
of being survivable in hostile environments (laser attack, 
particle impact), and of course, affordable to fabricate. Be- 
cause of the variety of extremely stringent, often unique, 
combination of material requirements, the only solution for 
many of these applications are novel advanced composites, 
which can be tailored to provide the required properties for 
specific SDI applications. 

Although, composite materials with high operating 
temperature potential have been under development for a 
number of years, there are major deficiencies in processing 
science and fundamental understanding of their mechanical 
behavior and environmental degradation mechanisms. 

In the past, because of their great complexity, composites 
have largely been developed in an empirical manner. As a 
result of the great combination of fibers, matrices, and com- 
posite architectures possible, as well as the near-term need for 


these materials, such inefficient and expensive material 
development procedures can no longer be tolerated. With this 
in mind, the Strategic Defense Initiative Office of Innovative 
Science and Technology established a program to establish 
scientific understanding of the processing science, composite 
behavior mechanisms, and new measurement techniques 
necessary to make possible the innovative materials required 
for future SDI systems. The program was designed to address 
processing on a fundamental level, identify and correct lack of 
understanding in specific areas, and demonstrate innovative 
material concepts. A center for processing science of metal 
matrix composites and ceramic matrix composites, two of the 
most important classes of materials for SDI applications was 
established during fiscal year 1985 at the Massachusetts In- 
stitute of Technology. In addition, other investigations were 
initiated to generate basic understanding of composite science, 
develop novel property measurement methods, or demonstrate 
innovative materials, at university and industrial research 
laboratories. Descriptions of the center and some of the re- 
search being conducted in carbon reinforced glasses, inter- 
metallic matrix composites, metal matrix composites, and on 
the development of new mechanical property measurement 
methods are presented herein. 
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A Center For The Processing 
And Evaluation Of Metal And 
Ceramic Matrix Composites, 


Materials Processing Center, 
Massachusetts Institute of Technology 


The great advantage of composite materials is the 
specificity of their properties, i.e., it is possible to engineer 
materials to a specific application. This advantage, however, 
has been a retardant to their application. With such freedom of 
design, it is not possible to accumulate a large data base. We 
may never have the production volume or the demand to 
produce such a data base for many systems. In the absence of 
a large data base, alternative methods must be found for 
predicting materials performance. Thus, the goal of the Center 
for the Processing and Evaluation of Metal and Ceramic 
Matrix Composites is the complete and fundamental under- 
standing of the relationships between microstructures, proper 
ties and processing parameters of metal and ceramic matrix 
composites. With this understanding, these technologies can 
be economical as well as a technical reality. 

The research topics sponsored by the SDIO/IST within 
the Center address issues pertaining to the successful introduc- 
tion of metal and ceramic matrix composites into engineering 
structures. Individual research topics have been selected on 
the basis of their ultimate impact on important SDIO structural 
applications and in their conformance with the program’s 
overall structure. 





Figure 1 
Specific Properties of A1 and Mg Matrix Composites 
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Metal and ceramic matrix composites have been simul- 
taneously addressed throughout the program. The two areas 
differ slightly in their principal foci: metal matrix research 
emphasizes improvements of property levels; while ceramic 
matrix composite research is concentrated on improved 
reliability. Concerns for the matrix-reinforcement interface, 
processing and microstructural control are fundamental to 
both. Although each activity is unique in its approach and 
goals, certain activities and approaches are common to all. The 





Central problem of metal matrix composites. Interfaces act as mechanical fuses to arrest impinging cracks by interface 
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emphasis of the Center is on fundamental process under- 
standing (i.e., understanding the relationships between 
processing, microstructure and materials properties). 

. Because of the low viscosity and short curing time (i.e., 
solidification time) of liquid metals as compared to resins and 
the net shape potential of casting technology, metal matrix 
composites are potentially less expensive to manufacture and 
more applicable to SDIO structural applications than resin matrix 
composites. That potential is being explored in the activities of 
the Center and the other programs associated with the Center. 

The initial work has been with Al and Mg alloy matrices 
because of the structural efficiencies of composites of those 
materials as shown in Figure 1. 

Carbon fiber (P-55) is being featured for most of the 
continuously reinforced studies. However, A1203 fiber in the 
form of Saffill™ and Fiber FP™ and SiC(AVCO SCS-2™) 
has been used in modeling studies. In the semi-solid slurry 
processing studies, SiC,B4C and A120s particulate materials have 
been used for either process development or process modeling. 

Carbon fiber reinforced composites are difficult to fabri- 
cate in any but the most simple shapes by any process other 
than casting or liquid metal infiltration processing. Further, the 
minimum sections required for the gossamer structures that are 
needed for space structures can not be produced by any other 
method. However, implementation of casting technology has 
been hindered by a lack of understanding of processing fun- 
damentals. These areas of needed knowledge include the 
following: 

Fluidity (fluid flow through a preform with or without 
accompanying solidification. 

macrosegregation. 

fiber wetting and interfacial reactions. 

the development of alloy microstructure during 
solidification fiber distribution 





Figure 3 


Matrix Transfer Molded Composite Pipe-Section 











Figure 4 


Composite Mirrors fabricated from Carbon Fiber Reinforced 
Glass 








Continuously reinforced metal matrix composites repre- 
sent the maximum potential for strengthening in directions of 
principal stresses. However, as much as a two-fold improve- 
ment in specific stiffness can be achieved isotopically with 
particulate reinforced metal matrix composites. The most 
economical method of producing discontinuously reinforced 
metal matrix composites is by solidification processing or, 
more specifically, semi-solid processing. As in pressure in- 
filtration mentioned above, the performance has been less than 
potential because of processing defects. The principal defects 
are porosity introduced during the inclusion of particulates 
into the melt to form the semi-solid slurry and poor distribution 
of ceramic reinforcement arising from the movement or “push- 
ing” action of the solidifying metal. 


— The critical fundamental issues that must be addressed 
are: 
The development of solidification microstructure in 
the presence of particulates. 
interfacial reactions and wettability. 
particle pushing effects. 
modeling of stirring and mixing 
modeling of deformation behavior of semi-solid slur- 
ries 





These issues are being addressed both fundamentally and 
practically. Fundamental studies are being conducted on the 
evolution of microstructure of Al-Cu-Mg and Mg-Al alloys 
with SiC and B,C particulates respectively. Particle pushing 
effects are being explored in Al-Cu-Mg alloy/SiC systems. 
Mathematical modeling of semi-solid stirring and rheology of 
deformation studies are also underway. An extensive cross 
cutting study is also underway to explore the wettability in 
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Figure 5 


Longitudinal and Transverse Strain Versus Applied Tensile 
load for 0-43 Vio HMU/7740 
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continuous and discontinuous systems as well as methods for 
treating fiber or particulate surfaces to improve wettability. 

The properties of composites are determined by the nature 
of the interfacial region. The need to engineer the interface is 
common to both metal and ceramic matrix composite ac- 
tivities in the Center. There are at least three reasons for 
tailoring (engineering) the interfacial zone in metal matrix 
composites; 1) to control the fracture path, 2) to protect the 
fiber surface from degradation due to handling or from chemi- 
cal interaction with the matrix during consolidation, and 3) to 
provide wettability for the ease of processing. Reasons 1 and 
2 are general for metal and ceramic matrix composites; reason 
3 is specific to solidification processing of metal matrix com- 
posites. 

The Center’s program is focused on tailoring the interface 
through use of Plasma Assisted Chemical Vapor Deposition 
(PACVD) This technique results in amorphous coatings with 
freedom to specify both coating modulus and residual stresses. 
The tailoring effort is backed by micromechanical analysis 
and measurement of specific properties of the interface as well 
as finite element elastic plastic analysis of fiber arrays. 

The central problem of engineering interfaces in metal 
(and ceramic) matrix composites is the arrest of cracks at the 
interface. In one program at the Center a P-55 carbon 
fiber/pure Al composite has been chosen as a model system 
for engineering the interface. The problem of a reinforcing 
fiber coated with a brittle sacrificial layer in a ductile matrix 
is shown in Figure 2. A crack is shown to nucleate in the 
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reaction zone between the matrix and the coating and impinges 
on the fiber/matrix interface. The stress Gee at6=0 probes the 
strength of the fiber while the same stress at6@="2 is responsible 
for interface separation. The radial shear stress 6,» at the 
interface (6=%2), on the other hand, probes the shear strength 
of the interface. Thus, two failure processes are possible upon 
interaction between impinging cracks and interfaces: 1) fiber 
fractures, or 2) fiber uncouples from the coating in shear or in 
tension. 

Interfacial strength requirements for the protection of 
fibers by crack arrest at the coating/fiber interface are given in 
Figure 2. When the ratio of; interface cohesive strength to fiber 
strength is less than the ratio of evolving elastic crack tip 
stresses responsible for interface separation and fiber fracture 
(i.e., the ratio of see atO=P2 see at@=0), the interface will 
separate in tension providing that the interfacial stresses are 
tensile. This condition leads to an upper bound for the interface 
cohesive strength. Likewise, when the ratio of interface shear 
strength to fiber strength is less than the ratio of Spe the 
interface (@=%2) to seg atO=0, the interface will separate in 
shear. This bounds the maximum interface shear strength. 
When the ratio of interface cohesive strength to interface shear 
strength is less than the ratio 5%,» at the interface, tensile 
separation will be preferred. Also, the toughness of the inter- 
face should be less than toughness of the fiber or the coating 
for the crack to advance through the fiber. 

Ceramic matrix composites have potential for a wide 
range of applications including components for advanced 
aerospace propulsion systems to laser hardening of space 
structures and orbital insertion and re-entry systems. Proper- 
ties generally sought include strength, hardness, stiffness, 





Figure 6 


Effect of Vol% Fiber on Composite Stress-Strain Behavior 
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Figure 7 


Tensile Fracture and Proportional Limit Strains for 
0?-HMU/7740 @ = Filament Tows, O= 3000 Filament Tows 
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Figure 8 


Fiber-Matrix Interfacial Shear Strength (tmax) vs Composite 
Short Beam Shear Strength (SBSS) for 0?-HMU/7740 with 
Various Matrix Additives 
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toughness, resistance to chemical attack, wear resistance, and 
reliability. Many applications require higher temperature 
capabilities than can be satisfied with alternative materials. 
The composite properties result directly from the properties of 
component materials, their bonding and the matrix/reinforce- 
ment architecture. 

The initial investigations in the ceramic matrix composite 
section of the Center have been directed toward reaction 
bonded SiC(RBSC) and reaction bonded Si3N«4 glass transfer 
molding (glass infiltration), and pre-ceramic polymer syn- 
thesis of BsC for fibers and matrices. Work is also in process 
on the evaluation of BN, non-stoichiometric silicon carbide 
(SiCi.x) deposited by PACVD on SiC fibers as well as the 
evaluation of the SCS-6 interface in SisN4 and SIC systems. 


Carbon Fiber Reinforced Glass 
For Space Based Applications. 


The use of carbon fibers to reinforce glass matrices 
provides a unique opportunity to achieve composites that 
possess properties particularly well suited to space based 
applications. As in the cases of reinforced polymers and me- 
tals, the use of carbon fibers provides composites with the 
potential for high levels of stiffness, strength, and dimensional 


stability. The use of glass as a matrix material offers an 
opportunity to achieve superior long time durability and resis- 
tance to environmental effects. Glasses are not likely to outgas, 
are highly resistant to atomic oxygen, and are exceptionally 
resistant to elevated temperature exposure. 

While a wide variety of ceramic matrices may also be 
considered as composite matrices, the use of glass provides the 
unique ability to fabricate composite articles of complex shape 
in a manner similar to that currently used for conventional 
thermoplastic matrix composites. An example of this 
capability is shown in Figure 3 where a thin walled, tube of 
woven carbon fibers in glass is shown. The tube was densified 
in a single operation with the actual matrix densification taking 
place in a matter of a few minutes. The resultant composite 
microstructure is fully dense. Significantly more complex 
shaped articles can also be fabricated as exemplified by the 
mirrors shown in Figure 4. These mirrors were fabricated with 
thin walled face and back plates encasing a lightweighted 
structural core. All elements of these mirrors were fabricated 
at United Technologies Optical Systems Division using carbon 
fiber reinforced glass. Their surfaces have been coated to 
provide the necessary reflectivity and they exhibit exception- 
ally high levels of dimensional stability. 

While the above describes several important advantages 
for the carbon fiber reinforced composite system, it can also 
be said that the true performance potential for such a system 
has not yet been fully realized. It is only the first of a broad 
range of fiber reinforced ceramic matrix composites which can 
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Figure 9 


Electron Microprobe Step Scan Analysis of Composites with 
Matrix Additives 
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revolutionize the performance of space based structures. The 
SDIO sponsored research activity is intended to create a 
fundamental basis for the understanding and improvement of 
composite performance and to lead to concepts for composite 
enhancement, not only for carbon reinforced glass but also for 
other future forms of ceramic matrix composites. 

The first issue addressed by this program is the unique 
Stress-strain behavior of ceramic matrix composites. As 
illustrated in Figure 5 for unidirectional carbon fiber rein- 
forced giass, the stress strain response to tensile loading is 
highly nonlinear. In addition, composite tensile strength is 
only 60% of that typical of polymer matrix composites. The 
nonlinear aspect of composite behavior is due to the 
microcracking of the glass matrix at a low value of strain 
and may prove to be the limiting factor in composite 
mechanical performance. The theory of Aveston, Cooper, 
and Kelly provides a basis for the determination of the level 
of strain at which microcracking is expected to take place. 
As shown in the following expression, increasing the 
volume percent of fiber in the composite should be effective 
in increasing this strain and as an additional benefit it 
should also increase composite stiffness. 


4 (1) 
— 12 T Ym Ef Ve |, 
Ec Eft Vm 


matrix failure strain in the composite 
elastic moduli of composite, matrix, and fiber. 
interfacial shear strength. 
= volume fractions of fiber and matrix. 
= fiber radius 
matrix surface fracture energy 
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Creation of composites with increasing levels of fiber 
content have confirmed this dependence of microcracking 
strain on fiber content, Figures 6 and 7. The strain level was 
increased five fold with an increase in fiber content from 25 
to 55%. 

As can be seen in the above expression, the fiber-matrix 
interfacial shear strength is also an important factor in control- 
ling microcracking. Increasing this level should enhance the 
microcracking strain and more importantly, it should increase 
composite shear strength as well. To demonstrate this, a series 
of unique glass compositions was created and used to react 
with the carbon fibers during composite fabrication. By these 
reactions it was possible to alter the nature of the fiber-matrix 
interface as confirmed by individual fiber-matrix debonding 
experiments. These have been related to composite shear 
strength with the result shown in Figure 8. As projected by the 
model, composite shear strength increased with increasing 
interfacial bond strength. Detailed chemical analysis of the 
interface region confirmed that reaction between the fiber and 
matrix (particularly in the case of the glass containing MoOs) 
was the cause of this enhancement. Electron microprobe 
analysis indicated that the Mo had migrated to the interface 
region (Figure 9), while x-ray diffraction showed that sig- 
nificant amount of Mo2C were created. 

The tailoring of the fiber-matrix interface properties af- 
fects many aspects of glass matrix composite performance. In 
addition to increasing composite strength, it is also important 
to preserve the interface as a region for resistance to crack 
growth in these brittle fiber-brittle matrix systems. At the 
present level of strength, it is possible to demonstrate tolerance 
to impact damage similar to that typical of resin matrix com- 
posites. As shown in Figure 10, only local damage is incurred 
by projectile impact. This is unlike the case of monolithic 
ceramics where total structural failure (shattering) would 
occur. 





Figure 10 


Front Impacted Faces of Graphite Fiber Reinforced Com- 
posites 














Figure 11 


Microstructure of NiA1 Composites (a) 0% TiB2 (b) 5% TiB2 (c) 10% TiBz2 (d) 20% TiB2 
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The current program is continuing to explore the fun- 
damental mechanisms controlling ceramic composite perfor- 
mance so that systems exhibiting superior high temperature 
strength and stiffness may also be created to meet the needs in 
space based power, reentry, and other applications. 


Nickel Aluminide/Titanium 
Diboride Composites 


by XD™ Synthesis 

Metal Matrix Composites are prime candidates for space 
applications in which high specific mechanical properties, 
high thermal and electric conductively, high vibration damp- 
ing and, dimensional stability are important. 

A new approach to forming metal matrix composites 
reinforced with XD™ particles not possible with standard 
metallurgical techniques, the XD™ process, has been under 
development at Martin Marietta Laboratories for the past 
several years. Under this program, an attempt is being made 
to apply the XD process to produce composites based on nickel 





Figure 12 


Transmission Electron Micrograph Showing TiB2 Particles in 
NiA1 Matrix 
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Figure 13 aluminides with enhanced high-temperature strength and low- 
temperature ductility. 

Room Temperature Hardness vs TiB2 Volume Content Of the various nickel aluminides, NiAl has many attrac- 
tive features, e.g., a very high melting point (1912 K)’, high 
modulus (higher than that of steel)®, low density (5.9 g/cm*)* 
and excellent oxidation resistance **, Its major drawbacks are 
low strength at high temperatures and inadequate low- 
temperature ductility.” Although grain refinement may partial- 
ly alleviate the ambient ductility problem® lack of adequate 
creep strength continues to be a disadvantage. The objective 
is to use XD™ synthesis to produce a dispersion of fine, 
discrete particles of a refractory second phase in this matrix, 
which will result in improved creep strength and enhance- 
ment? in toughness. An alternate approach to enhancing the 
high-temperature strength of NiA1 is to precipitate ternary 
compounds by adding group V elements.'° However, because 
these dispersoids are produced by precipitation, they will be 
unstable at high temperatures. The present approach will not 
suffer from this limitation since the dispersoid can be 
preselected to be thermally stable and has the added advantage 
that the matrix as well as the dispersoid can be formed in a 
single step. 

In the present investigation, NiAl with controlled volume 
fractions of titanium diboride ranging from 0.025 to as high as 
0.5 was made. Titanium diboride is of interest because of its 
high modulus, low weight, and good thermal stability. Sub- 
sequent x-ray diffraction (XRD) and scanning electron 
. “ microscopy (SEM) examination indicated that the thermal 
200; 10 15 stability of TiB2 in NiAl was excellent at least up to a tempera- 

ture of 1500°C. 
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Figure 14 
Flow Strength vs Strain Rate at 1200K 
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Figure 15 


Flow Strength vs Strain Rate 1300K 
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Figure 11 shows the microstructures of composites 
with 0, 5, 10, and 20 vol% tiB2; small amounts of A1203 
were present in all the samples. The average size of the TiB2 
particles was about 1 micrometer. Average grain size of the 
NiAl matrix decreased from about 30 micrometers in the 0 
vol% material to about 8-10 micrometers in the 10 vol% 
material which should result in enhanced toughness. Fur- 
ther grain refinement for second-phase contents exceeding 
10 vol% was found to be insignificant. A bright-field trans- 
mission electron micrograph from a composite with 5 
vol %TiBz is shown in Figure 12. The TiB2 particles appear 
to be well bonded to the matrix and show considerable 
substructure in the form of small angle boundaries, stacking 
faults and dislocation networks. The dark particles are 
A1203 measured densities for the various compositions 
agree well with calculated values, indicating adequate con- 
trol on the amount of second phase. Room-temperature 
hardness measurements made on composites with 2.5, 5.0, 
7.5, 10, and 20 vol% TiB2 (Figure 13) show a systematic 
increase in room-temperature hardness with increasing 
TiB2 content as expected. 

The high-temperature deformation behavior of these 
composites was evaluated by measuring their stress-strain 
response in air at 1200 and 1300K and strain rates ranging 
from 2.0 x 107 to 2.0 x 10“s". Their oxidation resistance 
was not adversely affected by the presence of TiB2. Flow 
strength vs strain rate data for composites with 0, 2.5, 7.5, 
10, and 20 vol% TiB2 are shown in Figures 14 and 15. The 
flow strength of composites with 20 vol% TiBz2 is about 
three times higher than that of the base composition. All the 
materials exhibited the expected dependence of flow 


strength on strain rate. The strain-rate dependence of flow 
strength for the composites is less than that of the base 
material. 

As mentioned earlier, a second important factor that limits 
the applicability of NiAl-based materials is lack of adequate 
low-temperature ductility. As tightly bound discontinuous and 
continuous composites have notoriously low fracture tough- 
ness, it is important to devise a manner to enhance the tough- 





Figure 16 


Damping Bahavior of Single-Ply (0°) P55/6061Al at Low Strain 
Amplitudes Indicating Little Strain Amplitude Dependence 
and Lower Damping than the 6061 Aluminum Matrix Alloy 
(Free-Free Flexure in Vacuum). 
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ness in these composites. One possibility would be to use 
transformation toughening"! of the NiAl matrix. A marten- 
sitic transformation occurs in NiAl at Ni-rich compositions 
and can be stress induced.’ If this transformation can be 
made to occur at the tip of a crack, toughness could be 
improved. Room-temperature toughness might also be im- 
proved by incorporation of weakly-bonded whiskers in the 
NiAl matrix. Both approaches are being pursued. It is an- 
ticipated that future efforts in the development of these 
composites will be based on the use of TiB2 for enhance- 
ment of high-temperature strength and the use of either 
transformation toughening or whisker-reinforcement for 
ambient toughening. 


Damping Characteristics Of 
Metal Matrix Composites. 


Primary concerns during acquisition, pointing and 
tracking (APT) maneuvers of large space structures (LSS) 
are to suppress vibration and achieve dynamic dimensional 
precision. Vibrational control can be achieved by various 
active control and passive damping measures. While the 
role of inherent material damping is frequently omitted in 
the active control analyses, it may provide a significant 
source of passive energy dissipation to ensure successful 
stabilization of large flexible space structures with fewer 
active controls.!* Inherent material damping capacity of 
conventional structural alloys is very low. Preliminary 


damping measurements of structural materials such as con- 
tinuous fiber-reinforced graphite-aluminum (Gr/Al) and 
graphite-magnesium (Gr/Mg) composites, which are impor- 
tant candidates for space applications, indicated that these 
metal matrix composites (MMC) exhibit high damping than 
their matrix alloys.'* In the current SDIO/IST program, a 
systematic investigation is being conducted to understand 
fundamentally the damping mechanisms and subsequently 
develop an approach to enhance the damping capacity of 
Gr/Al and Gr/Mg composites. A unique and innovative 
approach using in situ characterization techniques is util- 
ized to evaluate damping and to identify the sources of 
energy dissipation. 

Initially, reliable test techniques were developed to obtain 
repeatable and reproducible data of unidirectional single- 
ply P55/6061 Al (Gr/Al) composites. Three different test 
techniques, free-free flexure, clamped-free flexure (can- 
tilever beam) and tension-tension fatigue, were used. All 
flexural tests were conducted in vacuum because measured 
damping behavior can be severely affected by energy los- 
ses, particularly due to aerodynamic drag.'° Both flexural 
and extensional mode tests were conducted at strain 
amplitudes from 1 x 10° to 800 x 10° within the frequency 
range 0.1 to 150 Hz. Like metallic materials, Gr/Al com- 
posites exhibited regions of strain amplitude independent 
and dependent damping. At low strain amplitudes levels 
(e< 10°),the specific damping capacity, y, of Gr/Al com- 
posites was nearly independent of dynamic strain, and 
lower than the damping values of the aluminum matrix 





Figure 17 


Damping Behavior of Single-Ply (0°) P55/6061Al at Inter- 
mediate Strain Amplitude Dependent Behavior and Higher 
Damping that the 6061 Aluminum Matrix Alloy (Tension-Ten- 
sion Fatigue). 





Figure 18 


Comparison of Measured Damping Behavior versus the Be- 
havior Predicted by the Micromechanical Model (J2 Yield 
Criterion) for Single-Ply (0°) P55/6061Al Panels (Tension Ten- 
sion Fatigue). 
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Figure 19 


Granato-Lucke (G-L) Plot of Damping Data for Single-Ply (0°) 
P55/606 1Al under Tension-Tension Fatigue Testing Indicating 
the Proper Strain Amplitude Dependence. 
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Figure 16. Damping behavior in this region can be explained 
using Hashintheory'* which expresses the dynamic response 
to linear anelastic (or viscoelastic) composites in terms of their 
constituent fiber and matrix materials as shown below: 


wi = Vr (E/E )ve + (1-V_) (Em) Et/E11)Wm_— (2) 


where Wi1, Wr and Wm are the axial specific damping 
capacity of the composite, fiber and matrix respectively, Vr is 
the fiber volume and E};, Er and Em are the axial elastic moduli 
of the composite, fiber and matrix respectively. Generally, the 
damping contribution from Gr fiber is considered negligible 
as they are almost perfectly elastic. Using wr of about 0.1% 
and Wm of 1.2%, equation 2 yields yi: of 0.41%, which is quite 
consistent with the measured value of 0.6% by the free-free 
flexure method. Energy dissipation at these low strain levels 
can be attributed to anelastic relaxation mechanisms involving 
atomic diffusion, dislocation relaxation (Bordoni Peak), dis- 
location resonance and the Zener thermoelastic effect. 

At intermediate strain levels, the damping capacity in- 
creased with an increase in strain amplitude. Furthermore, the 
average for the composite was higher than the y value of the 
aluminum matrix, particularly in extensional mode Figure 17. 
Strain amplitude dependent damping response of these com- 
posites was examined in terms of both analytical 
micromechanical and microstructural models. 

A micromechanical model was developed (based on the 
assumption that the energy dissipation is caused by local 
yielding of the matrix at intermediate strain levels) using the 
commonly known “J2” yield criterion, i.e. effective damping 
energy depends upon distortional energy. This “J2” damping 


is anonlinear phenomenon and is affected not only by existing 
Stress (strain) levels, but also by stress and thermal histories. 
It was observed that at any dynamic stress level, the predicted 
roman y values correlated reasonably well with the measured 
damping values (Figure 18), provided that contributions were 
included from residual stresses (present near the interfaces due 
to differential thermal contraction between fiber and matrix 
during consolidation). These results indicate that the residual 
stress state plays a significant role in the damping behavior of 
Gr/Al composites. 

At the microstructural level, the strain amplitude depend- 
ent damping response was analyzed in terms of the Granato- 
Lucke (G-L) model.'”"* According to G-L theory, dislocations 
breakaway from minor pinning beyond a critical strain 
amplitude, and damping arises from the resultant hysteresis in 
the dynamic applied stress vs. dislocation strain curve. The 
strain amplitude dependent damping component wy can be 
expressed as w= (Ci/E) exp (C2/E) where C; and C2 are 
material constants related to microstructural parameters. Fig- 
ure 19 show a G-L plot of Gr/Al test data, and the calculated 
mobile dislocation density from this analysis is of the order 
10'? 1/m?. Dense dislocation networks near the fiber-matrix 
interface in Gr/Al composites were also revealed during TEM 
examination that can be attributed to the differential coeffi- 
cient of thermal expansion (CTE) between the fiber and 
matrix and the attendant residual stress state (Figure 20). 
Measured average dislocation densities of the order 10'¢ 
1/m?, obtained from various TEM images of inter-fiber 
matrix regions, is consistent with the mobile dislocation 
densities (10'* 1/m?) predicted by G-L theory, suggesting 
that dislocation damping is the operative mechanism at 
intermediate strain levels. 





TEM Image of Inter-Fiber Regions of Single-Ply (0° P55/606 1Al 
Indicating Dense, Tangled Dislocation Substructures Near the 
Fiber-Matrix Interface. 








Three/1989 31 





Damping test data analysis of Gr/Al composites showed 
that the matrix and attendant dislocation substructure near the 
fiber-matrix interfaces play a dominant role in energy dissipa- 
tion mechanisms. Among metallic materials suitable for com- 
posite fabrication, cast megnesium alloys such as Mg-1%Zr, 
Mg-1%Mn and Mg-1%Si exhibit very high damping in the 
range w ~ 50% - 70%.'? Subsequently, specimens fabricated 
from cast panels of P55/Mg-1%Mn and P55/Mg-1%2Zr are 
being tested to examine their damping response. Preliminary 
test results indicate that these Gr/Mg composites exhibit 
higher damping than the aforementioned Gr/Al composites 
(Figure 21). From these damping results and microstructural 
examinations, it is believed that damping can be further en- 
hanced by appropriate heat treatment to obtain the desired 
grain size and microstructure. 


Submicron Moire Interferometry: 
A Revolutionary Tool To 
Investigate The Micromechanical 
Response Of New Generation 


Materials 


Inherent to the concept of tailoring materials for a specific 
purpose is the ability to quantitatively understand the effects 
of such microstructural tailoring. The effects of tailoring a 
material spans orders of magnitude in effect - from neighbor- 
ing local constituents such as fibers and interfaces at the 





Figure 21 


Damping Behavior of As-Cast (0 °) P55/Mg-1%Zr Exhibiting 
Higher Damping than P55/606 1Al Composites (Clamped-Free 
Flexure in Vacuum). 
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micron level to gross structural level deformation at the cen- 
timeter level. This program is investigating the extension of 
an existing tool, Moire interferometry, into a revolutionary 
technique with an ability to visualize displacement fields over 
orders of magnitude yielding quantitative previously unseen 
phenomena from microstructural tailoring. The moire effect is 
observed when two superimposed repetitive patterns distort 
relative to one another and the distortion appears as a third 
pattern in space - a moire effect. Qualitatively, one terms the 
moire effect as parallel optical data processing. The moire 
effect is parallel due to the full field nature of the fringe pattern 
- it is not discrete. The data that is processed is the relative 
displacement between two gratings. The parallel nature of the 
moire effect is extremely powerful. If a grating is attached to 
a material such that it deforms with the material and is ob- 
served through an undeformed grating in space, the regions 
where material deformation becomes intense spatially be- 
comes evident in the full field moire pattern. The grating 
attached to the material is the active grating and it samples the 
material surface deformation. The grating that is held constant 
in space is termed the reference grating and it computes the 
deformation field from the sampling, active grating. The com- 
putation speed is the speed of light. Quantitatively, each moire 
fringe is a line of constant displacement whose direction is 
perpendicular to the direction of the reference grating. The 
value of displacement associated with a moire fringe is simple: 








U=N/F (3) 


where U is the displacement value of the moire fringe, 
N is the relative fringe order, typically from a datum. 
F is the reference grating computational frequency 
(lines/inch, 1/in, or lines/mm, 1/mm) 


It is evident then that the ability to sample material 
deformation is inversely proportional to the active grating 
frequency. The ability in accuracy of the computation is inver- 
sely proportional to the reference grating frequency. The upper 
limit of the mechanical moire line crossing effect is 1000 lines 
per inch (40 1/mm). This upper limit implies that the minimum 
deformation that could be detected on a material by a reference 
grating is 1/1000 of an inch (25 microns). More importantly 
the minimum spaual sampling that the active grating could 
sample is about 1/1000 of an inch. For new generation 
materials such as composites, the microstructural tailoring 
occurs on the micron level. Further the composites are ex- 
tremely stiff. Important deformations occur in micron size 
regions with typically a few microns in total deformation over 
an entire sample. The mechanical moire effect is inadequate 
for this class of composites. 

Moire interferometry eliminates both the low spatial sam- 
pling ability of the mechanical active grating and the low 
computational accuracy of the reference grating. In moire 





interferometry typical active grating frequencies are 30,000 
lines/in (1200 1/mm) with computational reference grating 
frequencies of 60,000 1/in (2400 1/mm). Moire patterns are 
being produced that can sample spatial detail in materials as 
low as 1/1200 of a mm (0.8 microns) with a sensitivity to 
compute in-plane deformation as low as 1/2400 of a mm or 0.4 
microns! By employing pulsed light sources we now have the 
ability to time sample transient deformation events as low as 
a millionth of a second (1 microsecond). What does this mean 
in physical terms? We now have the ability to look at the 
deformation fields from strong transient impacts in composite 
interfaces where high degrees of microstructural tailoring take 
place and where the subtle beginnings of material failure 
begin. 

Moire interferometry achieves such high grating frequen- 
cies by using the interference that occurs between two 
coherent beams of laser light interacting at a relative angle. 
The two beam interference yields a grating - whose frequency 


can achieve half that of the frequency of light employed. For 
blue - green light this upper limit is about 100,000 1/in or 4000 
1/mm! In moire interferometry the spatial grating is 
photographed and replicated onto the material. When the 
material is deformed in front of the two laser beam spatial 
grating moire patterns of extremely high sensitivity are 
produced. In addition, pulsed lasers may be employed to strobe 
light thus freezing transient deformation fields in advanced 
composites under impact. The laser pulse width is 20 billionths 
of a second (20 x10” sec.) - fast enough to freeze stress wave 
displacement fields in composite interfaces traveling in excess 
of 8 km/s! We are now beginning to understand the fundamen- 
tal mechanisms of transient impact delamination in com- 
posites. 

In Figure 22 moire interferometry was utilized to study 
microplasticity occurring between bundles of 10 micron 
graphite fibers and the magnesium matrix of a metal matrix 
composite subjected to vibration in a three-point bending 





Figure 22 


Static Moire Interferometry — Single Cycle Hystersis Damage. 
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mode. The discontinuity appearing in the displacement field 
is the magnesium microplasticity occurring at the fiber bundle 
interface - a mechanism of material damping previously un- 
seen. It is seen that the fringes really contain fine articulations. 
These fine articulations are individual 10 micron fibers shear- 
ing in the magnesium rich region that joins She larger bundles 
of 10 micron fibers. 

Figure 23 details the dynamic aspects of moire inter- 
ferometry. A graphite/epoxy composite is impacted to intro- 
duce transient waves parallel in the direction of the composite 
layering. This composite is composed of alternating layers of 
graphite fibers 90 degrees (perpendicular) to the impact direc- 
tion and 0 degrees (parallel) to the impact direction. The Figure 
shows fringes of constant displacement spanning the specimen 
free edge. Each fringe is an increment of compressive dis- 
placement of 0.4 microns. Within each spanning fringe are fine 
articulations corresponding to the number of alternating 0 and 
90 degrees layers. These articulations are due to the different 
speed of wave propagation between the 0 and 90 degree layers. 
Within the 90 degree layers the impact wave travels at 3 km/s. 
Within the 0 degree layers the impact wave travels at 8 km/s. 
The result is that these two different speeds meet at the 
interface causing a strong interfacial shear. Under stronger 
impacts delamination and spalling would occur. 


Concluding Comments 


Many of the future Strategic Defense Initiative systems 
are materials limited. These materials needs can best be met 
by tailor-making advanced composites, with properties op- 
timized for specific applications. 

The goal of the advanced composites program of the 
Innovative Science and Technology Office of SDI is to make 
possible new design options for design engineers developing 
systems for the Strategic Defense Initiative. The program 
approach is to demonstrate the feasibility of innovative com- 
posites and to establish a firm basis of scientific understanding 
of composite phenomena. 





Figure 23 


Dynamic Moire Interfermerometry — G/E Composite Stress Wave Behavior. 
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Abstract 


An existing technology is finding a new application in 
the study and development of laboratory x-ray lasers. 
Capacitively discharged pulsed power generators are being 
used to drive material loads to produce hot dense plasmas 
that possibly can be used as a medium for x-ray lasing or as 
a radiation pump source to photopump a variety of specially 
designed x-ray targets. The possibilities are considered of 
using pulsed power generators to drive a configuration of 
coaxially cylindrically symmetric krypton gas puff plasmas 
to conditions conducive to achieving gain in several lasing 
transitions of neonlike krypton. Using a fully self-consis- 
tent radiation hydrodynamics model, numerical simulations 
are performed driving the outer plasma inward where it 
crashes onto the core plasma establishing highly localized 
regions in the compressed core plasma. It is shown that in 
the localized regions it is possible to generate plasma con- 
ditions that can support a population inversion leading to 
gain in excess of unity in several transitions. Contour plots 
are presented during the collision process illustrating the 
behavior of the temperature, density, and gain in the core 
plasma. 
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Introduction 


Since the development and utilization of optical lasers 
there has been a steady flow of ideas and progress toward the 
development of lasers with shorter and shorter wavelengths. 
In the last few years this effort has been rewarded with the 
successful demonstration of lasing in the soft x-ray region of 
the spectrum. This accomplishment was achieved by irradiat- 
ing a specially designed target with a high powered infrared 
laser to create a hot dense plasma with properties conducive 
to lasing. The need to create a dense hot plasma as the material 
medium in which lasing will occur is due to the requirement 
that x-ray line emission is generally the result of atomic 
transitions between widely separated energy levels of highly 
charged ions. This suggests that the material medium cannot 
be in the solid or even in a neutral gaseous state but must 
instead be a hot plasma. This, of course, is the more conven- 
tional and time honored approach to x-ray lasers. It does not 
preclude some of the recent and exciting research involving 
free electron lasers, electron channeling, and multiphoton 
absorption concepts. However, these innovative approaches 
are just in their formative stages and need time to mature and 
reach their full potential. 





absorption concepts. However, these innovative approaches 
are just in their formative stages and need time to mature and 
reach their full potential. 

Since the majority of present day lasing schemes depend 
on the creation of a hot dense plasma, there are two separate 
but related issues that must be addressed. The first issue deals 
with the choice of a lasing scheme and a determination of the 
plasma conditions required to produce a population inversion 
and gain in the lasing transitions. The second issue deals with 
the selection of a power source capable of generating a plasma 
with the desired conditions. Both issues can be related by 
designing targets that are tailored to a particular power source 
in the hopes of producing a plasma with the desired conditions 
for lasing. Unfortunately, this is often easier said than done and 
generally requires a vigorous interaction and fine tuning be- 
tween theory and experiment. 

In general, the atomic processes leading to a population 
inversion and lasing in the x-ray region are the same ones 
responsible for lasing in optical lasers, only here they involve 
more energetic electrons in a hotter environment. Lasing 
generally occurs as the result of the microscopic processes 
between atoms (meant to include ions, also) and subatomic 
particles such as electrons and photons. The electrons are 
initially liberated from the atoms as a result of externally 
supplying substantial amounts of energy to the system. As the 
constituents of the system absorb the energy with a cor- 
responding increase in temperature, electrons are freed 


(ionized) from the atoms. The electrons now gain energy from 
the source and in their random walk through the medium crash 
into the more massive slower moving atoms and ions and 
further ionize the medium. The process avalanches and con- 
tinues until a kind of quasi-equilibrium is achieved and a 
balance is reached between energy in and energy out. During 
this process the average charge state of the plasma is increased 
and the plasma emits radiation; electrons are being liberated 
(ionization) and captured (recombination) and photons are 
being born (spontaneous decay and recombination) and either 
die (absorption) or escape the plasma carrying information 
about the environment it came from (in much the same fashion 
as astronomers learn about the stars). During this cavalcade of 
events the plasma percolates and its constituent atoms and ions 
have their levels populated, i.e., excited by the various col- 
lisional and radiative processes. A sample energy level or 
Grotian Diagram (as it is called) is shown in Figure 1 for a 
heliumlike atom (i.e., with two electrons orbiting the nucleus). 
Aselected number of levels are shown with their spectroscopic 
labels along with some of the processes involved in populating 
them. The process describing the evolution of level population 
is commonly referred to as ionization dynamics. The first step 
toward theoretically developing a laser is the requirement of 
creating a population inversion in one or more of the levels. 
An inversion occurs when the level population “piles up” over 
and above that normally achieved in most plasmas. The 
strategy to accomplishing an inversion is to choose an atom 





Figure 1 


Typical energy level diagram of heliumlike ion illustrating several processes leading to line emission. 
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Figure 2 
Schematic of Double Gas Puff 


Key energy levels of neonlike krypton. 
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whose atomic level structure, under the right plasma condi- 
tions, is conducive to overpopulation with respect to the 
ground state. 

We will focus on the neonlike system of krypton (i.e., 
krypton in an ionization state with only ten bound electrons 
remaining) because its unique atomic structure makes it a 
prime candidate for establishing a population inversion in the 
manifold of excited states. Once the inversion is obtained, then 
the competition between the various atomic processes 
depopulating the inverted states occurs. If radiative decay to 
selected lower states dominate the de-excitation process, the 
emission of a monochromatic beam of soft x-rays will follow. 
Gain in the lasing transitions will occur if the intensity of the 
X-Fays in transit through the material is described simply by 
the relationship I = Ip exp (GI) where I is the emitted intensity 
of the beam at a distance 1 for its origin with intensity I,.The 
gain coefficient, G, is positive and depends on the atomic and 
plasma properties for the medium. Under ordinary circumstan- 
ces, a beam of light traversing a medium of length 1 is 
attenuated along its path length and emerges reduced in mag- 
nitude, depending on the opaqueness of the medium, i.e., G is 
negative and represents the absorption coefficient. For gain to 
occur the wave must be exponentially amplified rather than 
reduced in magnitude while traversing the medium. In prin- 
ciple, with the appropriate physics models and numerical 
simulation codes, it is possible to theoretically determine 
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under what conditions inversion and lasing will occur and 
whether such conditions can be created in the laboratory. We 
now turn to the issue of generating the plasma. 

The power source considered here has been around for 
some time now and initially was employed extensively in the 
early days of magnetic fusion research, viz., pulsed power. 
More recently, pulsed power drivers have been used success- 
fully to generate plasmas for use as intense high brightness 
sources of x-rays for nuclear weapons effects testing. Concep- 
tually, the simplest of these devices is the so-called z-pinch. 
The z-pinch in its basic configuration is a cylindrically sym- 
metric plasma carrying an axial current which both heat and 
confine the plasma by its own induced magnetic field. As the 
current increases, the magnetic field also increases and com- 
presses the plasma, accelerating it radially inward, causing 
increased heating by converting the kinetic energy of inward 
motion to thermal energy; containment and heating occur 
together. Unfortunately, when driven by very large currents 
these plasmas exhibit pathological behavior by becoming 
unstable and ultimately breakup. For both fusion and x-ray 
laser applications this behavior is obviously undesirable. On 
the other hand, one of the most attractive features of the 
z-pinch plasma is its natural elongated shape as a medium in 





which to exhibit lasing and gain. Whether it is possible to use 
these plasmas in an x-ray laser scenario is actively being 
pursued. 

It recenily has been suggested that by properly design- 
ing material loads it should be possible to delay or deter the 
growth rate of the instabilities for times longer than the 
characteristic lasing time.” Preliminary experimental 
results based on concentric gas puffs appear to support this 
hypothesis.’ In order to take advantage of these new find- 
ings on stability our target is configured as concentric 
axially symmetric krypton gas puffs. The outer gas puff is 
a hollow cylindrical annulus while the inner gas puff is 
uniform and symmetric about the axis. The scheme is con- 
ceptually shown on Figure 2. Discharging a megajoule 
capacitor bank through the gas puffs causes rapid ionization 
by ohmic heating. As the current continues to flow and 
increases to a peak value of between 3 to 5 magampers the 
induced magnetic field stops the plasmas’ outward expan- 
sion and begins to compress the plasma accelerating it 
radially inward on a collision course with the inner plasma. 
The dynamics of the implosion is often characterized as a 
“slowplow” because it behaves in a similar fashion by 
scooping up and compresssing more and more mass as it 
moves along. During the run-in the outer plasma can 





Figure 4 


Schematic of Numerical Simulation Code. 
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achieve speeds in excess of a million miles per hour before it 
crashes and stagnates on the inner plasma. During this phase 
oftheimplosiontheinductiveormagneticenergy isconstantly 
being converted to kinetic, thermal, ion internal, and radiative 
energy. 

Because the plasma evolves through a series of 
temperatures and densities, its radiation pattern similarly 
evolves from the infrared to the visible, the ultraviolet, and 
the x-ray region of the spectrum, becoming a broadband 
radiator. In order to produce a plasma with just the right 
temperature and density history for lasing to occur, it is 
important to match the drivers elecirical characteristics 
with the load design and configuration. Since we are focuss- 
ing attention on the neonlike ionization stage the electrical 
parameters must be chosen so as to optimize the plasma 
conditions for its production. Once into the neonlike stage 
the level populations are determined by the ionization 
dynamics. An abbreviated Grotian energy level diagram is 
shown in Figure 3 along with the level’s spectroscopic 
identification, quantum numbers, and the wavelengths of 
the lasing transition in Angstrom units (10°'° meter). In 
addition, the spontaneous decay rates of the lasing transi- 
tions are also shown on the Figure; the notation, e.g., of the 
172 A line between the 3p and 3s states is shown as 2.26'° 
which translates to 2.26 x 10'° sec."! giving the lifetime of 
the upper 3p state as 1/2.26 x10"° sec."! = 4.42 x 10" sec. 





Initial conditions of Gas Puffs. 
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Figure 6 


X-ray emission spectra including line and continum com- 
ponents. 
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Model 


The implosion dynamics of the krypton gas puff plasma 
is described by the one dimensional equations of radiation 
hydrodynamics. For simplicity we ignore the magnetic field 
and replace its effect on the system by assuming the outer 
plasma has an initial kinetic energy and moves with a 
constant velocity. This replicates the initial stages of the 
implosion dynamics but avoids the complex chain of events 
that occur during the first nanoseconds when the initial 
surge of current encounters the gas puffs causing phase 
changes and eventually ionizes them. Also, since there is 
reasonable evidence from the preliminary analysis of recent 
experiments and numerical simulations that lasing, if it is 
to occur, will be achieved well before the onset of the 
destructive interference caused by the formation of plasma 
instabilities, we will ignore the problem of plasma in- 
stabilities. With these caveats, the model contains enough 
of the relevant physics to provide a reasonable repre- 
sentation of events describing the behavior of the plasma 
column as a lasing medium. 

The overall morphology of the evolving gas puff plasmas 
can be schematically represented by appealing to the adjacent 
flow chart, Figure 4, characterizing the system of equations 
used to describe the chain of events. 


40 Naval Research Reviews 


In the diagram the different boxes denote the various 
physical processes involved in characterizing the implosion 
dynamics. The arrows entering and leaving each box indicate 
physical quantities serving as input or output data for each of 
the subroutines. Although the structure of the flow chart is 
suggestive of a scaler process, the actual calculation is non- 
linear and requires a series of interactive iterations between 
the various subroutines before a unique solution is converged 
upon. The actual form of the equations, numerical techniques, 
and numerical zones are presented elsewhere (for those inter- 
ested in the details of the simulation see reference 4). The 
equations of hydrodynamics, denoted as HYDRO, describes 
the fluid properties of the plasma including shock wave 
phenomena. TCON and EOS are ancillary subroutines repre- 
senting thermal conduction processes and equation of state 
information, respectively. The lower portion of the schematic 
diagram represents the ionization and radiation physics. In- 
cluded are the ionization dynamics of nonequilibrium plas- 
mas, radiation absorption and emission processes, and the 
transport of radiation. All these processes, from above and 
below, are funneled into a kind of “mixing bowl” labeled 
Interface. The mixing bowl assimilates all the information and 
passes the quantities up and down and iterates on them until a 
converged set of quantities is obtained in each zone of the 
numerical mesh. 

A number of computer simulations were carried out for a 
variety of gas puff conditions. One of the more interesting 
simulations resulted from the initial set of conditions shown 
on Figure 5. Both inner and outer puffs are composed entirely 
of krypton gas with about seven times more mass contained in 
the outer puff (140 gms compared with 20 gms). The 
distribution of mass is Gaussian with the bulk of mass located 
at 15 mm and 1 mm, respectively from the origin. The initial 
temperature of the outer (inner) gas puff plasma was taken to 
be 2.Sev (0.Sev). A tenuous background plasma pervades the 
space between the inner and outer plasmas and has a density 
of 5x10” gm/cm?. This background level was chosen for 
numerical convenience and is sufficiently low that it has no 
influence on the results. The outer plasma starts moving 
inward with a speed reaching upwards of about 1 million 
miles per hour (actually 4.5 x 10’ cm/sec.). The leading 
edge crashes into and stagnates onto the inner plasma. The 
“collision” continues in a snowplow-like fashion with more 
and more mass being accumulated and stacked up on the 
inner plasma. During this collision process the kinetic ener- 
gy of run-in is redistributed into thermal, ionization, radia- 
tion and kinetic energy associated with the generation and 
propagation of shock waves and the outward motion of 
material after the “bounce.” This oscillatory motion con- 
tinues until it is damped by vicious forces. The interface 
region between interacting plasmas is where most of the 
interesting phenomenon occurs and where we expect lasing to 
be observed. 








Figure 7 


Plasma Temperature as a function of time and radius. 
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Figure 8 


Plasma lon Density as a function of time and radius. 
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Figure 9 


Fraction of neonlike ions as a funtion of time and radius. 
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Figure 10 


Gain per cm of the 170A line as a function of time and radius. 
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Figure 11 


Gain per cm of the 172A line as a function of time and radius. 
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The radiation emitted during this sequence of events not 
only cools the plasma but is the messenger of conditions within 
the plasma. The emission spectra is a “fingerprint” that can 
lead to knowledge about the temperature, density, size, 
opacity, and velocity of the emitting region. A typical emission 
spectra is shown in Figure 6 at 36.33 nsec into the implosion. 
The snapshot picture of the radiation signature shows the 
intensity of radiation in terms of spectral lines (spike-like 
structures) sitting on top of a continuum of radiation composed 
of recombination and bremsstrahlung as a function of photon 
energy extending from 10 - 10,000 eV. Most of the line 
radiation from 400 to 2,000 eV originates from neonlike 
krypton. This is very encouraging because it is in the neonlike 
system that we expect lasing to occur. The temperature as- 
sociated with the spectra is several hundred eV and the particle 
density of neonlike ions is between 10'* and 10” ions/cm? in 
the interface region. The curves represented by T=637 and 247 
correspond to a plasma radiating as a Blackbody at those 
temperatures in eV and are simply there to indicate what a 
blackbody spectrum would look like in comparison with the 
actual radiated spectrum. 

A novel way of presenting the results of the numerical 
simulations is to illustrate the behavior of the physical 
parameters such as temperature, density, and gain as a function 
of time and radius in the region where the two plasmas interact. 
In the following Figures time runs along the vertical axis from 
27 nsec to 36 nsec and displacement or radius runs along the 


horizontal axis from the origin out to 0.08 cm. In this region 
the core plasma is being compressed and heated by the incom- 
ing outer plasma. Initially, the core plasma is relatively cold 
but as the outer plasma gets hot it starts radiating with some 
fraction of the photons impinging on the core plasma causing 
its temperature to rise. This process continues until the bulk of 
the outer plasma stagnates onto the core plasma causing fur- 
ther heating. During the collision the core plasma becomes hot 
and begins to expand, radiate, cool, and recombine. The 
various shades of light and dark are coded on the Right Hand 
Side (RHS) of the figures. There are several ways to under- 
stand and interpret the results. One of the easiest is to choose 
a time and then observe how the parameters vary from larger 
to smaller displacements and simply repeat the process for 
increasing times. Also, remember that these results extend into 
the plane of the figure for the length of the cylinder, which is 
3 cm, or for a distance over which the quantities are uniform 
and homogeneous, which appears to be about 1.7 cm?>, The 
radial temperature history is shown on Figure 7 over distances 
mainly confined to the compressed core plasma and the inter- 
action region. The brighter regions are primarily the result of 
compressional wave heating or weak shocks. The expanding 
darker regions later in time extending from the axis outward 
reflect the effects of radiative cooling while the brighter 
regions from 34 nsec on, indicate heating due to weak shocks 
moving outwards. The temperature peaks around 1 kilovolt 
(over 11 million degrees) in this region. The ion density is 
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shown on Figure 8 where it is seen that as time increases the 
density along the axis increases where the plasma is compress- 
ing. These high pressures near the origin slow down the 
incoming plasma creating velocity gradients at about 33 nsec. 
Thus, kinetic energy is constantly being thermalized in this 
region, commencing just after the on-axis assembly. Radiation 
cooling in this region falls off sharply with radius following 
the corresponding density falloff. Just after 33 nsec, the plasma 
near the origin begins to expand in a complicated fashion; parts 
of the plasma move outwards while other parts continue their 
inward motion. This produces compressions and rarefactions 
leading to heating and cooling in localized regions. At just 
after 34 nsec, the plasma at about 0.008 cm from the origin is 
expanding radially while elsewhere the plasma is moving 
radially inwards. From the contour plot of temperature shown 
on Figure 7, it can be seen that alternating rapid cooling and 
heating is produced from 33 to about 36 nsec in localized 
regions. Very large pressure gradients which have been built 
up at 24 nsec will very quickly be relieved. At 35 nsec relief 
of the large pressure gradient has caused a virtual stagnation 
to occur at about 0.035 cm. Plasma is now moving radially 
outward from this point and colliding with imploding plasma. 
As a result, a fairly broad shell of hot plasma (about 1 kev) has 
been produced. Since the ion density is fairly low over much 
of this region the radiative cooling rate is not excessive. 

The fraction of neonlike krypton ions in the ground state 
is shown on Figure 9. Since the fractional abundance of 


neonlike ions is strongly temperature dependent and only 
weakly density dependent, it is not to surprising to observe that 
the number of neonlike ions tracks the temperature. The 
greater abundances occur for temperatures in excess of about 
400 ev which is about the threshold temperature for the ap- 
pearance of neonlike krypton. Shown on Figure 9 are anumber 
of localized regions where the fractional abundance of neon- 
like ions exceed 50% and a few regions where it is upwards 
of about 70%, i.e., in those regions of the plasma anywhere 
from 50 to 70 percent of all the ions are in the neonlike ground 
state. A very propitious situation, indeed, for our proposed 
x-ray laser scheme! 

By way of an introduction to the final gain calculations a 
few words are in order explaining the relationship between the 
shape of a spectral line and the gain in the lasing transitions. 
Whenever a transition between discrete states of an atom ina 
plasma takes place a spectral line is radiated from the plasma. 
Contrary to the results presented on Figure 6 spectral lines 
emitted from plasmas do not appear as infinitely sharp lines 
or spikes. Under closer scrutiny with good wavelength resolu- 
tion it is possible to observe a spectral line profile as a function 
of wavelength is generally distributed symmetrically about its 
unperturbed wavelength representing the separation of the two 
atomic energy levels involved in the transition. The line profile 
or distribution function results from such line broadening 
effects as the finite lifetime of the upper state of the atom, 
which is a result of the Uncertainty Principle, collisions and 





Gain per cm of the 184 A line as a function of time and radius. 
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environmental influences which also manifest themselves in 
broadening the spectral line. Oftentimes it is just these effects 
on spectral lines that provide valuable information on condi- 
tions within the environment. In the gain calculations 
presented here the line profile functions are represented by 
Voigt functions, i.e., a convolution of Lorentzian and Doppler 
width profiles where the damping parameters for the 
Lorentzian is composed of natural and collision broadened 
widths and the Doppler includes thermal and motional 
broadening effects. The transport of radiation in the presence 
of large velocity gradients is accounted for in the Sobolev 
approximation developed for the expanding envelopes of 
Stars. 

Conditions conducive to a population inversion leading 
to gain in the lasing transitions were achieved for three of the 
lines: the 172, 170, and 184A lines, respectively. These results 
are shown on Figures 10, 11, and 12, respectively. Reference 
to Figure 3 shows the wavelengths of the respective transi- 
tions. Gains of one or more occur in highly localized regions 
eventually achieving values in excess of 9/cm for the 172A 
line; 5/cm for the 170A line; and 3/cm for the 184A. The 
atomic processes contributing to the gain are predominantly 
due to recombination and electron impact excitation. The gain 
172A line extends over the largest radial region extending 
from the origin out to about 0.048 cm from 28 to 29 nsec. A 
more localized region occurs from 34 to 36 nsec. The gain for 
the 170A line occurs at the same time but is confined to smaller 
regions. Finally, gain for the 184A line, which is due 
predominantly to electron collisions, is confined to an even 
smaller region peaking around 35 nsec. Note, that on the scale 
on the RHS of the gain profiles some of the values represented 
are negative numbers and correspond to absorption of radia- 
tion. The point to be made is that it is possible to achieve a 
population inversion and gain in three of the expected lasing 
lines from our imploding plasma scenario. The final confirma- 
tion of whether this approach is feasible is currently unfolding 
with the analysis of the experimental observations. Prelimi- 
nary results are very encouraging. 
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Artifact Diamond 
Its Allure and Significance 


Max N. Yoder 
Office of Naval Research 


Abstract 


While the preponderance of the mechanical, optical, and 
electronic properties of natural diamond have been known for 
over a decade, only recently has artifact diamond in tech- 
nologically useful form factors become an exciting possibility. 
The advent of sacrificial, lattice matched crystalline substrates 
provides the basis not only for semiconducting applications of 
diamond, but for optical mirrors, lenses, and windows as well. 

As a semiconductor, diamond has the highest resistivity, 
the highest saturated electron velocity, the highest thermal 
conductivity, the lowest dielectric constant, the highest 
dielectric strength, the greatest hardness, the largest bandgap 
and the smallest lattice constant of any material. It also has 
electron and hole mobilities greater than those of silicon. Its 
figure of merit as a microwave power amplifier is unexcelled 
and exceeds that of silicon by a multiplier of 8200. For 
integrated circuit potential, its thermal conductivity, saturated 
velocity, and dielectric constant also place it in the premier 
position (32 times that of silicon, 46 times that of Gas). 
Although not verified, its radiation hardness should also be 
unmatched. 

Aside from its brilliant sparkle as a gemstone, there has 
been little use of diamond in the field of optics. Processing of 
the diamond surface now appears to be as simple as that of any 
other material — albeit with different techniques. In fact, it 
may be possible to etch diamond far more controllably (at 
economically viable rates) than any other material as the 
product of the etch is gaseous and the etched trough is self- 
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cleaning. Other properties of diamond make it an ideal optical 
material. Among them are its unmatched thermal conductivity, 
its extremely low absorption loss above 228 nanometers, and 
unmatched Young’s modulus, Poisson’s ratio, tensile strength, 
hardness, thermal shock, 2nd modulus of elasticity. If the 
recently-found mechanisms by which erbium impurities in 
III-V junctions can be made to “lase” are indeed applicable to 
indirect gap semiconductors, then even efficient diamond 
lasers at attractive portions of the spectrum become a distinct 
possibility. 


Background 


Although diamond has been known for centuries, it was 
the advent of X-rays that permitted its crystalline structure to 
become known in the second decade of the present century. 
The father and son Bragg duo used X-ray diffraction to show 
that diamond had a cubic structure and a lattice constant of 
3.56 Angstroms. With tetrahedral bonding, the nearest neigh- 
bor distance is 1.53 Angstroms and at the centerpoint the 
electron density is approximately 1.59 electrons per square 
Angstrom! Although the covalently bonded diamond crystal 
exhibits a very open lattice, the presence of an interstitial 
impurity atom, a misplaced carbon atom, or a vacancy leads 
to a considerable change in the electron density and to lattice 
distortion. 





Characterization 


With an electronic bandgap of 5.45 electron volts, 
diamond should exhibit no optical absorptions within the 
visible spectrum as its absorption edge should be 2250 
Angstroms. Although many water clear gemstones are found 
in nature, there are also those of yellow, blue, green, pink, 
mauve, and black. Clearly nature has placed electronic absorp- 
tion levels within the bandgap. It was these in-gap absorption 
levels that led diamonds to be broadly classified as type 1 or 
type 2. If the specimens demonstrated a sharp increase in 
absorption at wavelengths shorter than 2250 Angstroms and 
there was negligible absorption at wavelengths longer than 
3000 Angstroms, they were classified as being type 2. Type 2 
diamonds exhibit very little visible fluorescence under 
moderate ultra violet or ionizing illumination. Under intense 





Figure 1 


Optical absorption in diamond and other materials. 
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ionizing radiation, even type 2 diamonds, will exhibit lumines- 
cence in the blue region. Although ultraviolet emission has 
been seen, it is very rare. In type 2 diamonds, the response to 
radiation near to but below the bandgap is generally that of its 
response to radiation of much longer wavelengths. In this case 
the dielectric constant throughout the visible spectrum is equal 
to the static dielectric constant (5.76) or the square of the index 
of refraction. Type 2 diamonds are further subdivided as A or 
B with the A group having a very sharp and temperature 
independent absorption cut off (they are P-type semiconduc- 
tors) and the B group having a broader and temperature 
dependent absorption cut off. Type 1 diamonds exhibit sig- 
nificant absorption at wavelengths of 3000 Angstroms or 
longer. Exception for their tribological use as abrasives or 
cutting tools, they have little commercial usefulness. Although 
diamond is monatomic and non-polar, even the best of natural 
diamonds exhibit some infrared absorption centered at 7.8 
micrometers. Even then, its absorption is approximately 4 
order of magnitude less than that of competing materials as 
shown in Figure 1.! 


Mirrors, Lenses, and Coatings 


Figure 1 compares the optical absorption of diamond with 
several competing materials for optical lenses and windows. 
Although not shown in the figure, diamond also exhibits 
approximately four orders of magnitude less peak absorption 
than do magnesium difluoride, yttrium oxide, and zinc 
selenide. For use with high intensity optical beams such as 
generated by free electron lasers, diamond has another un- 
matched attribute; its thermal conductivity is unexcelled. Thus 
even if appreciable energy is absorbed, the absorbed energy 
can be readily dissipated to the edge of the diamond. As such, 
optical beams exhibiting power densities greater than 10° 
Watts/cm? can be transmitted by diamond windows.” 

Liquid etches, even at high temperatures, are well known 
to be ineffective in etching diamond. While atomic oxygen is 
known to attack diamond, it is not an efficient etch. Using a 
combination of energetic ions and a physisorbed, oxygen-con- 
taining gas, Geis has shown that the diamond surface can be 
etched as easily as that of any other semiconductor.’ Not only 
has he been able to create etch troughs ten time deeper than 
the remaining diamond pillars are wide, the residual effluent 
is gaseous and the trough bottoms are self-cleaning! The 
feasibility of a diamond Fresnel lens is thus established. Using 
special high power laser processing techniques, Geis has also 
demonstrated that diamond windows can be fabricated as 
optical polarizers giving extinction ratios in excess of 15:1 at 
temperatures as high as 1400 Celsius.‘ 

While it is difficult to conceive of the use of natural 
diamonds to protect optical mirrors, lenses, and windows 
made of other materials, the advent of synthetic diamond 
engenders the optical engineer with better approaches to sys- 
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tems design. Not only can the underlying optical material be 
protected against microparticle abrasion by virtue of the un- 
matched hardness of diamond (9000 kG/mm?), but its ex- 
tremely high surface energy (as high as 9400 ergs/cm?) 
protects it from many (including biological) adherents — 
virtually nothing (excepting hydrogen) wets diamond. To im- 
prove the threshold of high power optical damage to silicon 
windows, a consortium of investigators from Crystallume, 
NASA Langley Research Center, and Old Dominion Univer- 
sity investigated synthetic diamond overcoatings.° They used 
the thermal stress parameter to motivate their work. It is 
proportional to tensile fracture strength, thermal conductivity, 
thermal coefficient of expansion, and elastic modulus — in 
each of which the properties of diamond are unexceeded. The 
measured threshold of single pulse 1.06 micrometer illuminat- 
ing damage ranged from a low of 5.97 J/cm? to a high of 12.43 
J/cm*. Comparisons of the thermal stress parameter of 
diamond with that of other optical materials is shown in 


Figure 2. 


Electronically Active Diamond 


Virtually all of the electronics properties of diamond are 
unexcelled excepting that of electron mobility. Even then, its 
mobility exceeds that of silicon. With very high charge carrier 
mobilities and velocities, diamond is expected to exhibit ex- 
tremely fast electronic devices capable of operating with very 
large voltage excursions, at very high temperatures, and with 
high immunity to radiation damage.° Unfortunately, diamond 
is an indirect bandgap semiconductor and it is well known that 
injection lasers fabricated with indirect bandgap materials 
have not been successful. Diamond, however, affords a new 
approach. With its extremely large forbidden bandgap, there 
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exists the possibility that many elements of the periodic table 
may have optical transition energies therein, as blue emission 
and even externally pumped lasing is well known in diamond.” 
In recent years, charge carrier pumping of impurities in more 
common semiconductors has produced emission and even 
single mode lasing.*° Using similar techniques it is feasible 
that diamond can be internally pumped by forward biasing a 
diamond junction or, more plausibly, a heterojunction of boron 
nitride on diamond. Since diamond is a column IV semicon- 
ductor, a valence 3 rare earth impurity could serve as the 
acceptor material immediately inside the P side of such a 
junction and possess a very large capture cross section of 
“pump” electrons. Figure 3a illustrates such a rare earth atoms 
shown as A3+ before forward bias while figure 3B illustrates 
a possible electron and optical transition after application of 
the forward bias. In this illustration 4f transitions of rare earth 
atoms are shown, but 2d transitions of the transition metals 
impurity atoms may also be considered. Figure 4 illustrates 
transition levels from impurity atoms having energy levels just 
within the conduction and valence bands. While less probably 
than those transition shown in Figure 3, they are possible. 
Finally, optical transitions from impurity levels just within the 
conduction band to energy levels within the same ion, but 
within the diamond bandgap can not be ignored. Thus it is 
conceivable that a diamond injection laser could be engineered 
to emit at virtually any wavelength from the infrared to the 
ultra violet. 


Summary 


The combination of this tensile strength, unexcelled ther- 
mal conductivity, low optical absorption, very large bandgap, 
appropriate Young’s modulus and Poisson’s ratio, and recent 
advances in its synthesis render it an ideal candidate for many 
optical applications. 





Figure 3 


Intra ion energy level transitions in diamond. (a) before forward 
bias, (b) after bias pulse 
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Bandgap magnitude intra ion energy transitions 
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The Oxygen Heat Engine 
Based On Solid 


Electrolytes 


A. V. Joshi, N. Weber, & J. Nachlas 
Ceramatec Inc. 
Salt Lake City, Utah 


Abstract 


The oxygen heat engine, based on solid state ceramic 
oxygen ion conductors like ZrO2 and Bi2Os, is a promising 
thermoelectric energy conversion concept being developed at 
Ceramatec, Inc. The device converts heat to electricity direct- 
ly, at high efficiency. The device utilizes thermal energy from 
sources such as solar concentrators, and expands oxygen gas 
across a ZrO2 ceramic membrane at 900°C or greater and 
compresses the gas through BixO3 membrane at about 500°C; 
electricity is generated electrochemically across the ceramic 
membranes. The actual efficiency under load is expected to 
exceed 20. Several uniqne features such as the absence of 
moving parts, low mass, compactness and modularity, high 
efficiency, and simplicity of the system make this concept very 
attractive for space platform and satellite applications. 
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In a preliminary research and development effort spon- 
sored by ONR/SDIO, the feasibility of this novel concept has 
been established and a bench scale heat engine has been 
successfully built and tested. For realistic commercial 
development of such a device, additional engineering is re- 
quired. Such an effort would involve further development in 
the area of electrode/electrolyte systems that are durable, 
efficient, and withstand the operating temperatures. Addition- 
ally, work towards integrating sub-systems to provide a multi 
mega-watt power system capable of delivering useful power, 
with the high integrity necessary for critical mission require- 
ments, is required. Such improvements have been made in 
recent technology developments of high power solid oxide 
fuel cells. 





This article describes the concept of the oxygen heat 
engine. Recently acquired experimental data that 
demonstrates “proof of concept” are described. The ad- 
vantages over present state-of-the-art technologies are pointed 
out. Future directions for research and development are dis- 
cussed. 


Introduction 


Current space vehicles and platforms are generally 
powered by solar photo-voltaic generation systems coupled 
with storage batteries. These systems offer adequate power for 
satellites requiring relatively low power levels; typically about 
10K We/- operation. Future space power applications, how- 
ever, demand much greater power due to different mission 
requirements. Tomorrow’s space platforms are envisioned to 
require 300 K We or greater by the mid to late 1990s; some will 
periodically consume many megawatt-hours of energy.’ Such 
space platforms will include manned space stations, com- 
munication stations, surveillance platforms and defensive 
weapons. Clearly future space power needs are several orders 


of magnitude greater than what can be effectively delivered by 
current technologies. The challenge is to explore new power 
technologies that will meet the projected power needs. A 
potential solution is the use of solar or nuclear power based 
heat engine systems. 

About twenty years ago, scientists at the Ford Motor 
Company Scientific Laboratory conceived the idea of produc- 
ing electric power from a thermally regenerative 
electrochemical system called the sodium heat engine (SHE), 
based on a sodium ion conducting solid electrolyte.” This 
invention generated a lot of interest in the scientific com- 
munity due to the significant advantage offered by solid 
electrolyte technology. The advantages, over current state-of- 
the-art technology, are no moving parts, modularity, compact- 
ness and light weight, and high system efficiency. 

In pursuit of these advantages, Ceramatec scientists con- 
ceived of a solid electrolyte energy conversion device based 
on an oxygen ion conducting solid electrolyte. This device, an 
oxygen heat engine (OHE), offers some advantages over the 
SHE. In the following sections, the working principle of the 
oxygen heat engine is explained. The projected performance 
of this device is compared to the SHE, and preliminary 





Figure 1 


Ceramatec'’s oxygen generators based on zirconia tubular design are shown in the photograph. Ceramatec demonstrated over 


5,000 hrs. of continuous high efficiency performance through these generators. Model A is 15 liters/hr. capacity unit, while Model 
B is 3 liters/hr. capacity unit. 
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laboratory performance data are described. Finally the ap- 
plicability of the oxygen heat engine as a potential power 
source for future space power applications is discussed. 


The Oxygen Heat Engine 


Basic Concept 


The critical materials that make the operation of this novel 
technology possible belong to a class of materials called solid 
electrolytes. These materials have the unique property of being 
able to transport ions rather than electrons through their crystal 
structure. Zirconia (ZrO2), and bismuth oxide (BizO3) are 
examples of oxygen ion conducting solid electrolytes. 

Although the concept of using a ZrO2 based solid 
electrolyte ceramic for pumping oxygen and for an oxygen 
heat engine has been around over a decade** the first practical 
small prototype, using single tube cells, was built and tested 
by Ceramatec in 1982 to develop and evaluate an efficient 
electrolyte/electrode system. The initial experiments in 1982 
at Ceramatec on ZrQ>-based tubular configurations 
demonstrated the concept but unfortunately excessive power 
losses were experienced due to electrode/-electrolyte deficien- 
cies. In 1984, Ceramatec developed a novel type of 





Schematic diagram of electromechanical heat engine based 
on solid oxygen ion conductors. 
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electrode/electrolyte system which resulted in an order of 
magnitude lower power consumption for the zirconia oxygen 
generator. With this development, the oxygen ion conducting 
solid electrolyte process became attractive in terms of com- 
mercially viable oxygen pumping devices. 

In 1985 Ceramatec built many small prototype tubular 
oxygen generators with capacities ranging between 1/2 liter 
and 10 liters/h of oxygen flow. These prototypes demonstrated 
over 5,000 hours of continuous, efficient performance under 
current densities exceeding 100 mA/cm”. A photograph show- 
ing these devices is given in Figure 1. 

In Ceramatec’s oxygen heat engine, gaseous oxygen is 
expanded at a high temperature across a ZrO2 membrane and 
compressed at a lower temperature across a BizO3 membrane, 
much like a conventional heai engine. The use of ZrO2 for 
electrochemical compression at low temperature is not ap- 
propriate because the conductivity of ZrO2 at 500-600°C is too 
low to achieve any useful power density. This fact was recog- 
nized by Ruka & Weissbart® who patented the use of a ZrO2 
membrane only for the high temperature expander in an 
oxygen heat engine. They did not recommend the use of aZrO2 
membrane for the low temperature compressor but instead 
suggested that the compressor might be “any of the well 
known mechanical devices for compressing gases”. It would 
appear that an all electrochemical device without any moving 
parts using oxygen as the working fluid has been undeveloped 
for want of an efficient oxygen ion conductive 
electrolyte/electrode to serve as an electrochemical compres- 
sor at temperatures low enough so that (T2-T}) is sufficiently 
large. The recent developments of high oxygen ion conduc- 
tivity in modified Bi2Os electrolyte materials.”"!* at tempera- 
tures as low as 500°C presented an opportunity to use these 
materials for an electrochemical compressor. Ceramatec 
scientists proposed this concept to SDIO and it is illustrated in 
Figure 2. As shown in Figure 2, oxygen gas is expanded 
through a ZrO2 membrane at perhaps 1000-1200°C (T1) and 
compressed through a modified BizO033 membrane at 500- 
600°C (T:). During operation, oxygen at a high pressure P2 is 
expanded electrochemically and isothermally to pressure P}. 
This EMF at zero current can be expressed in terms of Nernst 


(1) 


where R = gas constant 

F = Faraday constant 
Part of this output voltage from this cell is used to drive an 
electrochemically similar isothermal BizO3 based membrane 
compressor cell at a lower temperature. The net output voltage 
Vo is the difference, in the Nernst potentials at T2 and T; i.e. 


Vo = (RT2/4F) In (P2/P1) - (RT1/4F) In (P2/P1) 
= [R (T20-T201)/4Fy in (P2/P1) 
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Heat engine cycle comparison 
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The molar heat of expansion, q, abstracted from the 
reservoir at T2 is (RT2) InP2/P;. The theoretical maximum 
efficiency E is just the ratio of net work out (4 FV)/heat in at 
T2 per mole of fluid transferred. 


E = 4FVo/q = (T2-T1)/T2 (2) 


The cycle also includes the sensible heat absorbed and 
given up as oxygen circulates between the two thermal reser- 
voirs. If we imagine a perfectly efficient gas-gas heat ex- 
changer for the two fluid streams, then the two terms of equal 
magnitude and opposite sign cancel out in the expression for 
efficiency. 

As with any heat engine, the larger the T2-T1, the larger is 
the efficiency. The upper practical temperature limit for this 
oxygen based device is governed by electrode stability. The 
practical lower temperature is governed primarily by 
electrolyte conductivity, but also by other thermally activated 
processes associated with mass and charge transfer at inter- 
faces. In Figure 3, P-V and T-S diagrams are shown for the 
oxygen heat engine cycle. Shown also for comparison are 
corresponding diagrams for several other gas cycles. The 
similarity between the OHE cycle and the Ericsson cycle is 
evident. According to Karlekar'? a gas turbine plant with 
reheat, intercooling and regeneration executes a cycle that is 
close to the Ericsson cycle. 


How Does Oxygen Heat Engine 
Compare With Sodium Heat 
Engine? 


The oxygen thermoelectric generator has many attractive 
features in common with the sodium heat engine for SDI 
applications. One of the clearest potential advantages of this 
solid electrolyte based electrochemical device is that it can be 
very small and light weight. There are no moving parts except 
for slowly circulating fluid. Modular designs are possible 
allowing for system expansion without major redesign or 
scrapping of existing hardware. Smaller generators will be 
equally economically feasible. The cell is simple, rugged and 
may be insensitive to ionizing radiations. Efficiency will be 
virtually independent of size. The nature of this electrochemi- 
cal technology allows significant design flexibility to meet any 
specific end use requirements. 

Compared to the sodium heat engine (SHE) for example, 
the oxygen heat engine has some advantages. A heat exchange 
between input and output legs of the high temperature reser- 
voir is easier to accomplish. Because the working fluid is an 
insulator, achieving an electrical series system to minimize 
lead thermal conduction and resistive losses will be much 
easier. A gaseous fluid system will also be less sensitive to 
thermal shock and acceleration forces. Finally, the relative 
ease of working with oxygen over sodium may reflect in 
capital and maintenance cost advantages. The comparison also 
reveals some disadvantages. Two membranes are required in 
the cycle instead of one, and the oxygen solid electrolytes, 
even at high temperature, are inferior in conductivity to 
sodium beta"-alumina. The working fluid cannot serve as a 
conducting lead, and four permeable electrodes are required 
instead of one. The output voltage of the expander cell is 
smaller by a factor of four for a given pressure ratio since four 
electrons are required to transfer one mole of fluid compared 
to one in the case of sodium. 

An analysis of the oxygen fluid device will undoubtedly 
show maximum efficiencies at a lower electrolyte area for 
specific power levels than are found for the SHE. Available 
from previous studies on the SHE, the analytical and theoreti- 
cal machinery for a detailed calculation of thermal conversion 
efficiencies in the oxygen system is ready to be put in motion. 


Experimental Results 


The program at Ceramatec has focussed upon materials 
development, fabrication process development, device set up 
and assembly, and performance testing and evaluation. 

The major thrust of uur initial efforts was directed towards 
the understanding, synthesis and development of low tempera- 
ture BizO3 based membranes, as well as efficient high 
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Figure 4 


Oxygen ion conducting solid electrolyte tubes with various 
compositions. 
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temperature and low temperature electrode systems. While the 
high temperature electrolyte membrane has been widely 
studied, the low temperature bismuth oxide membrane needed 
development especially due to the thermodynamic instability 
exhibited by these electrolytes under actual operating condi- 
tions. 

As pointed out by several investigators”! Bi2033 has 
four polymorphs a, p, y, & 5 BizO3. The 5 phase exists be- 
tween 1000 K and 1097 K. This phase has an FCC defect 
fluorite structure and was found to have the highest oxygen 
ion conductivity with a low activation energy. Attempts to 
stabilize the 6 phase at lower temperatures, under non- 
electrochemical conditions, were successfully made by 
various scientists and are well documented in the literature.~"! 
Asuccessful oxygen heat engine device must not only contain 
a high conductivity phase Bi2O3 electrolyte membrane for the 
efficient compressing of O2 at low temperatures, but it must 
exhibit exceptional durability and phase stability in an 
oxygen/vacuum environment and under an oxygen pumping 
mode. We have carried out stability studies on these materials 
under actual pumping modes. We will now describe the 





Figure 5 


Schematic of bielement cell Bench Scale Apparatus 
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Stability Studies for 2 Electrolyte Systems Current Decay as a Function of Time Under a Constant Applied Potential of 0.2 Volts 
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fabrication of components, test apparatus as well as ex- 
perimental results obtained on a bench scale bi-element heat 
engine device. 


Fabrication. 


BizO3 Electrolytes: Bi2O3 electrolyte tubes stabilized 
with various aliovalent dopants such as Ba, Pb, Y, and Nb 
nominally 12mm dia, 25cm long and 1.1mm wall were fabri- 
cated at Ceramatec by cold isostatic pressing followed by high 
temperature sintering. The processing conditions were studied 
to achieve concentric, straight axis, low porosity (95 theoreti- 
cal density) gas impermeable tubes. These tubes are shown in 
Figure 4. The processing conditions were optimized to 
produce these tubes. 

ZrQ2 electrolytes: Stabilized zirconia tubes nominally 
25cm long and 1cm dia were fabricated by the usual procedure 
of cold isostatic pressing followed by sintering at high 
temperature of 1600°C for ? hours. 

Electrodes: For initial experiments, various noble metals 
as well as Ceramatec’s proprietary electrodes were used with 
Bi2Os electrolyte based on our prior extensive experience with 
electrochemical devices such as oxygen pumps. These 
materials are good electrodes for feasibility testing of the 
concept for operation up to 500 hours at temperatures below 


700°C. These electrodes were applied by painting the 
electrolyte tubes (both ZrO as well as Bi2Os) using a nylon 
brush, followed by baking at 700°C to 900°C to drive off the 
organic binders in the electrode slurries. 

Bench Scale Device: A simple tubular bench scale bi-ele- 
ment device (Figure 5) containing a ZrO2 ceramic expander 
and a Bi2O3 based electrochemical compressor was fabricated 
and used for feasibility demonstration- Parametric experi- 
ments varying the composition of the BizO33 compressor cell 
and its temperature, while keeping the high temperature ZrO2 
cell at a constant temperature of 900°C were conducted using 
this bench scale apparatus. Stability and device performance 
characteristics were studied and are reported below. 


Stability of Bi203 Electrolytes 

As pointed out previously, building a successful oxygen 
heat engine device depends critically on utilizing the high 
conductivity phase of BizOs electrolyte that must be stable 
under operating conditions. Stability tests under dynamic 
electrochemical conditions were carried out on Bi2O; tubes 
with the various dopants as outlined below. Experiments con- 
sisted of applying 0.2 volts D.C. on the respective electrolyte 
tube fitted with electronically conducting electrodes which 
were kept heated at 700°C. The current resulting from the 
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applied potential is directly related to amount of oxygen being 
pumped through the solid electrolyte membrane. As an ex- 
ample, 1 ampere of current corresponds to 3.5 cm?/min. of 
oxygen gas transport through the electrolyte. The current was 
monitored as a function of time over two hundred hours. Our 
experiments indicate that certain compositions of Bi2Os, al- 
though possessing high conductivity, were unstable under 
(oxygen pumping) testing conditions. For example, despite the 
fact that the ionic conductivity of BizOs(Y203) compositions 
is high at temperatures as low as 600°C our experiments 
indicate that these materials are not useful as solid electrolytes 
under the application of voltage in air environment. We have 
found that Y203 doped Bi2O3 composition decomposes into 
the lower conductivity rhombohedral and monoclinic phases. 
Furthermore, it was found that significant grain growth oc- 
curred in these solid electrolyte tubes under applied potential 
conditions. It appears from our experiments that the phase 
diagram of Bi203.Y20s as given by Datta and Mahan’* is 
insufficient to interpret the behavior under applied D.C. cur- 
rent conditions. Some doped Bi2O3 compositions (such as 
Group II element doped compositions) exhibit a rapid decay 
in current carrying capacity within 24 hours. No further inves- 
tigations were carried out on such compositions due to the 


severe degradation of these electrolytes. The aim of our 
detailed investigations was to identify those BizO3s composi- 
tions that underwent negligible degradation and hence led to 
stable device performance. To illustrate device performance 
decay, two compositions of doped Bi2O3 were investigated for 
current density measurements as a function of time under an 
applied potential of 0.2V. The current decay curves are shown 
in Figure 6. Preliminary analysis indicated that the current 
decay in these materials may have occurred due to degradation 
in electrode/electrolyte interface as well as electrolyte insta- 
bility. Further experiments are in progress to characterize the 
electrode/electrolyte interface and identify electrochemically 
stable Bi2Os electrolyte materials. Our recent laboratory find- 
ings show promise for the development of such stable 
materials. 


Bench Scale Test Results 


As described previously, the high temperature element 
consisted of a stabilized ZrO2 tube with electronically con- 
ducting electrodes, while the low temperature cell was a 
circular shape of stabilized Bi2O3 material. In order to simplify 
the analysis of the data, the high temperature cell was always 





Figure 7 
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kept at 900°C and the temperature of the low temperature cell 
was varied. Two compositions of Bi2O3 solid electrolyte were 
studied, namely BizO3(Y 203) and BizO3(Nb20s). 


Figure 7 shows the power v. load characteristics for Y203 
stabilized as well as Nb2Os stabilized electrolytes at two different 
temperatures. These data clearly indicate that BizO3.Nb20s solid 
electrolyte is superior to BizO3. Y 20s solid electrolyte in terms of 
power density. The voltage vs. current curves for both the 
electrolytes are shown in Figure 8. Further studies are underway 
to characterize other stabilized BizO3 based materials for use in 
the low temperature compressor in terms of their stability and 
influence on the power output of the heat engine. 


Future Developments 


The current program has clearly demonstrated the concept 
of oxygen based thermoelectric electrochemical conversion 
utilizing zirconia as a high temperature expander and Bi2O3; as 
a low temperature compressor. However, for achieving the 
high power densities, the performance characteristics have to 
be optimized. This will be achieved through improvements in 
understanding of the electrode and electrolyte materials as 
well as the fundamental charge transfer processes that occur 


at the electrode/electrolyte interface. Theoretical analysis of 
the heat engine with attention given to mass transport con- 
siderations will be conducted to further assess the performance 
characteristics. 

Further development can also be made via design and 
fabrication technology refinements. Virtually all rapid advan- 
ces occurring in the Solid Oxide Fuel Cells (SOFC) technol- 
ogy (Ceramatec’s Planar Geometry SOFC and Argonne 
National Laboratory’s Monolithic SOFC design) can be ap- 
plied to oxygen heat engine. With such thin film designs and 
efficient electrodes, power densities of over 1000 mW/cm? 
are possible. 

The program has also demonstrated that for commercial 
viability of this device, further effort in the development of 
efficient, durable & stable electrode/electrolyte systems is re- 
quired. For an electrochemical system based on ZrO? and BizO3 
solid electrolytes, this entails safe, reliable, and efficient operation 
for a period of time in excess of 40,000 hrs at operating tempera- 
tures of 500-700°C for BizOz and 1000-1200°C for ZrO2 based 
electrochemical cells under 1 atmosphere of oxygen pressure and 
lower. When a durable and efficient electrode/electrolyte system 
is demonstrated, advanced engineering will be required. Efforts 
will then focus on design and economic issues to allow more 
complete and realistic assessment of the oxide ion solid 
electrolyte concept for thermoelectric power generation. 





Voltage -vs- Current for various Low Temperature Solid Electrolyte Elements 
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With the capability and experimental results acquired in 
the present program, we believe that oxygen ion conducting 
solid electrolyte approach offers unique opportunities and 
exciting possibilities for space power generation systems. 
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Research Notes 


A New Optoelectronic Semiconductor Switch -— Op- 
toelectronic semiconductor switches have received consider- 
able attention, for use in high power switching. However a 
drawback of this device is that it requires a high power laser 
source operating for the entire duration of the current flow 
through the switch. A new concept being researched by Profs. 
K. Schoenbach and K. Lakdawala at Old Dominion Univer- 
sity, offers promise to overcome this problem. The material 
used for switching is GaAs doped with deep level impurities 
such as copper. The switching mechanism utilizes the deep 
levels as intermediate energy levels for charge storage. The 
increase in conductivity is accomplished by electron ioniza- 
tion from deep centers and bound hoies. The reduction in 
conductivity is obtained by hole ionization from excited 
centers and subsequent recombination of free electrons and 
holes. The entirely new and unique capability of this switch is 
that it uses photons to turn-on and turn-off the conductivity on 
command, but there is no radiation required to sustain it during 
the on-state. (Dr. Gabriel D. Roy, ONR) 


Radiation-Hardened Non-Volatile Semiconductor 
Memory - Transistors with Silicon-Oxide-Nitride-Oxide- 
Silicon (SONOS) gate structures that can act as the memory 
elements for a Non-Volatile Semiconductor Memory (NVSM) 
have been demonstrated by Lehigh’s Prof. Marvin White. The 
work, which was done under a contract entitled “Reliability of 
Smail Geometry VLSI Devices for Microelectronics”, is 
funded by SDIO and managed by ONR. The devices are 
rad-hard to at least 5 MRad total dose, and are capable of 
greater than 10’ Write/erase cycles, greater than 10'? read 
cycles, and retention time greater than 1 year. The trapping 
kinetics must still be studied and modeled, but the results 
obtained thus far suggest that near term development is war- 
ranted. The potential impact lies in the use of these devices in 
large scale, reliable, rad-hard, NVSM arrays that are needed 
for a variety of DOD applications. The first use of this tech- 
nology could be in the SDIO Boost Surveillance and Tracking 
System (BSTS) program. (Dr. A. Goodman, ONR) 


Direction-Of-Arrival of Signals Using Radon Trans- 
forms — Professor R. L. Kashyap, an ONR contractor at 
Purdue, supported through SDI funds, has proposed a new 
procedure to estimate the directions-of-arrival of emitters. The 
procedure helps to locate the emitters from signals received 
by an array of sensors in an electronically disturbed and noisy 
environment. The two dimensional data is transformed into a 
set of one dimensional sequences by the Radon Transform and 


processed to locate the peaks of the two dimensional spatial 
spectrum by a robust technique. The procedure does not re- 
quire any prior information about the source emitters and does 
not assume any model structure for the unknown colored 
noise. The algorithm is highly suitable for parallel processing. 
(Dr. R. N. Madan, ONR) 

Reliability Predictions of VLSI Circuits —- During a 
recent SDI/IST program review, Professor Hu of the Univer- 
sity of California at Berkeley presented his latest research 
results on VLSI reliability. Since metal electromigration is one 
of the primary causes of failures for next generation VLSI 
devices, he has developed a model to predict the lifetime 
reliability based on AC stress testing. The model predicted that 
AC stress lifetime would be 10,000 times longer than DC 
stress lifetime. This prediction has recently been verified by 
laboratory experiments. The model improves the accuracy of 
prediction on the lifetime of VLSI thin-line metals by several 
orders of magnitude over the present predictions. He also 
presented some results on charge trapping and oxide break- 
down in submicron VLSI devices, which are also reliability 
concerns. His research is a very important part of trying to 
improve the reliability of VLSI chips and systems, which are 
widely used in many Navy systems. (Dr. C. Lau, ONR) 


Radiation-Harden Nonvolatile Computer Memory - 
One of the requirements for space-based electronic systems is 
radiation-harden nonvolatile memories. Professor M. White 
of Lehigh, under ONR/SDI/IST funding, is developing such 
in reliable, small VLSI devices. The new device called 
SONOS (Silicon Oxide Nitride Oxide Semiconductor), with 
a gate dielectric consisting of very thin multilayers of silicon 
oxide, silicon nitride, and silicon oxide. These SONOS 
devices have been shown to be resistant to radiation of over 1 
Megrad, and can be used in nonvolatile memory applications 
with retention time well over 10 years. Together with radiation 
harden VHSIC technology, this SONOS technology when 
fully developed may prove to be extremely important to 
Navy’s space-based memory-intensive satellite systems. (Dr. 
C. Lau, ONR) 


New Compiler for Rapid Antenna Array Weight 
Design — Professor Steven Boyd, recipient of an ONR Young 
Investigator award (1985-88) at Stanford, has been investigat- 
ing linerization of nonlinear circuits and systems. He has 
recently applied his techniques to computer aided design of 
antenna array weights. Antenna array weights are variable 
signal amplitude and phase angles assigned to each array 
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element that give an improved array beam pattern. Variable 
antenna element weights improves antenna performance in 
terms of desired mainlobe width, maximum side-lobe con- 
Straints, position of nulls, and percentage of bandwidth 
coverage. 

His compiler translates these antenna specifications into 
a non-differentiable convex optimization problem and deter- 
mines the feasibility of the specifications. If feasible, a com- 
puter program calculates the optimal weights needed to 
provide the desired antenna pattern. Many examples have 
been worked for planar, circular and parabolic antenna array 
geometries of up to 30 elements. The results have been sub- 
stantially better than other computational methods. Dr. Boyd’s 
work has generated considerable interest in the Navy with a 
NWC, China Lake collaborative effort in progress. (Dr. R. 
Madan, ONR) 


ONR Scientific Officer is invited Speaker at IEE Con- 
ference in London — Dr. Gabriel D. Roy was an invited 
speaker at the Pulsed Power Conference sponsored by the 
Institution of Electrical Engineers, UK. The conference was 
held at the Savoy Place in London. Dr. Roy presented his paper 
on “High Power, High Repetition Rate Switches: An Over- 
view”, which was then published in the IEE Proceedings. 
While at UK, he visited the universities where US 
(SDIO/ONR) supported research programs are carried out. 
These were Aston University in Birmingham, Univ. of Kent 
in Canterbury, Univ. of Strathclyde in Glasgow and Univ. of 
St. Andrews in St. Andrews. The UK research focusses on 
surface breakdown phenomena in dielectrics and insulators, 
diagnostic development of surface breakdown and the impact 
of meteorites in space. He attended the annual program review 
of the SKI supported pulse power programs conducted at the 
ministry of Defence, London and was a special guest at 
Imperial College, London at the Executive Committee meet- 
ing of the UK Space Power Group. 


Diamond Growth Speeds Up — Research Triangle In- 
Stitute (RTI) with SDIO/ONR support, has demonstrated a 
method of growing diamond at very fast rates without the 
expenditure of excessive amounts of energy. The new method 
parallels one recently advanced by a group at AT&T Bell 
Laboratories and uses alternate cycle epitaxy wherein pure 
methane and molecular hydrogen are alternately applied to the 
growth surface — both in the presence of excited state helium. 
In this method it would appear that during the hydrogen 
portion of the cycle, the diamond copiously absorbs atomic 
hydrogen and effusively releases it to the growth surface 
during the methane portion of the cycle. In this manner there 
is effectively no boundary layer to inhibit the arrival of atomic 
hydrogen to the surface and growth rates of 50 to 100 times 
those of conventional low energy approaches are experienced. 

While this method is expected to be applicable to the 
growth of single crystalline diamond for electronics use, the 
preliminary studies have all been conducted on polycrystalline 
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material where it has a near term application in optics. Lenses, 
windows, and mirrors made of these films are expected to 
operate in the infra-rad, visible, and ultra violet portions of the 
spectrum where their optical losses are lower than those of 
competing materials and their ability to resist microparticle 
abrasion and scratching is unsurpassed. 

With a coefficient of friction similar to that of Teflon, 
extremely hard diamond films are expected to provide durable 
and friction-resistant bearings. As a semiconductor material, 
diamond is superior to silicon or gallium arsenide for high 
power electronics applications. Workers within this field an- 
ticipate that within six years diamond films will be available 
for use as semiconductors at a cost equivalent to present 
gallium arsenide materials. As a semiconductor, diamond is 
unexcelled in virtually all of the properties which characterize 
semiconductors. Only in electron mobility is diamond inferior 
to gallium arsenide; however, it remains superior to silicon. 
Mobility is represented by the slope of the lines in the figure. 
The properties illustrated translate to a figure of merit for 
microwave power amplifiers that is 900 times that of gallium 
arsenide and 8000 times that of silicon. The ability of semi- 
conducting diamond to operate at very high temperatures 
makes it an ideal candidate for use in monitoring jet engines 
and nuclear reactors. (Dr. Max Yoder, ONR) 
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